Method and arrangement for analyzing the interaction of high frequency electromagnetic emissions wit

专利摘要:
Method and arrangement for analyzing the interaction of electromagnetic high-frequency immissions (HFI) with vegetative regulatory mechanisms of a test subject, comprising the following method steps: determining the course of the power density of HFI over a defined period, determining the variance of the time interval by means of electrocardiogram measuring method (ECG ), consecutive heartbeats and the course thereof determined autonomic function parameters (HRV-FP), synchronization and calibration of the HFI and HRV-FP curves to a reference time interval (T1_n), compared with characteristic, HFI-induced secondary significances, which have been determined from the current experimental design or from at least one of these preceding experimental design with the current test subject or third test subjects and in the form of digital or analog data, preferably in the form of graphic gradients, alphanumeric or algorithm mix data stored on a storage device; upon detection of a defined, at least partially present agreement of the currently ascertained HRV-FP course with characteristic, HFI-induced secondary significances or on correlation with primary significances in the HFI course, a positive evaluation of the presence of its humanbiologically relevant influencing of the vegetative regulatory mechanisms of the test subject takes place ,
公开号:AT516204A1
申请号:T666/2014
申请日:2014-08-29
公开日:2016-03-15
发明作者:Michael Matissek；Peter R Mmag Hauschild
申请人:Peter R Mmag Hauschild；Michael Matissek；
IPC主号:

专利说明:

The invention relates to a method for analyzing the interaction of radio frequency electromagnetic emissions (HFI) with vegetative regulatory mechanisms of a test subject according to claims 1 and 2 as well as to corresponding arrangements for carrying out the method according to the invention according to claims 29 and 30.
Wireless information technologies using technically generated electromagnetic fields have been widely used in recent years. Electromagnetic fields (EMF) emitted by high frequency sources serve as carrier waves to which respective information signals are modulated. Common transmission standards are GSM, UMTS, LTE, WLAN, Bluetoooth / IEE802, DECT and others.
When an electromagnetic wave hits an absorbing body, the energy contained in it is converted to another state or heat. This thermal mode of action is also based on the limits set by the International Commission on Non-Ionizing Radiation Protection (ICNIRP 1998) for the protection of the general public from EMF exposure. Their recommendations were adopted by the EU (Council Recommendation 1999) and the ÖVE / ÖNORM E 8850.
In addition to the aforementioned thermal effects, athermal effects of radiofrequency electromagnetic fields, which correlate with adverse health effects, have also been found in a large number of scientific studies (Santini 2002,
Navarro et al. 2003, Hutter et al. 2006, Abdel-Rassoul et al., 2007, Blettner et al. 2008). From a study published in 2009 at the Institute of Environmental Hygiene at the Medical University of Vienna, it is estimated that adverse effects on the human biology, such as those caused by irradiation flux density of 0.5 mW / m2 (see also the ÖNORMadmissible thermal limit values of 4,500-10,000 mW / m2). Increased immunosuppression is causally detectable (Hutter, Kundi, 2009). As a precautionary measure, the Salzburg Medical Council even recommends an electromagnetic radiation flux density of 0.001 mW / m2.
Despite empirical evidence of athermal EMF effects, their mode of action with respect to the biological organism has not yet been conclusively clarified. For example, the formation of free radicals or oxidative and nitrosative stress is discussed as a plausible mechanism of action at the intra- and intercellular level (Friedmann et al., 2007, Simko 2007, Pall 2007, Bedard and Krause 2007, Pacher et al., 2007, Desai et al. ,
Since the development of high-frequency technology can only look back on a relatively recent history, there are no long-term studies in this area per se. Until such findings are obtained from such field studies, many regulatory and environmental practitioners recommend the greatest possible reduction in EMF exposure for precautionary reasons. For example, following an in-depth review of all currently available pertinent studies, the WHO International Cancer Research Agency IARC noted in June 2011 the need to reduce high frequency EMF and classed it as "possibly carcinogenic to humans" (Group 2B).
The precautionary principle has already been taken into account, for example. by the Swiss Federal Council with a "Non-ionizing Radiation Protection Regulation (NISEV)" for places of sensitive use ("OMEN") such as living quarters, schools, hospitals and retirement homes, with greatly reduced EMF plant limits compared to the regulatory regime. EMF-related complaints are recognized as physical restriction in some countries, such as the UK and Sweden, under the name Electrohypersensitivity (EHS), along with the statutory rights of EHS-affected workers to reduce exposure.
In the course of occupational health care, an increasing number of companies decide to voluntarily reduce electromagnetic pollution. For example, the Bavarian has
Automobilkonzem BMW restricted the permissible high-frequency electrosmog exposure at the workplaces of its approximately 105,000 employees worldwide to an in-house maximum of 100 pW / m2, commonly known as the "BMW limit value". In order to keep this limit, physical measures for exposure reduction have been taken in office and laboratory buildings, e.g. DECT and WLAN base stations shielded with metallised glass panes and wireless routers equipped with attenuators. Companies implementing such screening programs record a lower number of sick leave and more job satisfaction among their employees.
The finding of EMF-related complaints is difficult in practice and presents environmental physicians with extensive differential diagnostic requirements, as set out in the "Guideline for the Examination and Treatment of EMF-Related Disorders and Diseases" of the Austrian Medical Association of 3 March 2012. While individuals do not show immediate symptom symptoms due to their individual ability to regulate, other already relatively low levels of EMF exposure respond significantly. If EMF sensitivity is present, exposure reduction is indicated in each case.
Although previous studies on the influence of high-frequency EMF on EEG or ECG or heart rate values of a subject have revealed correlations of such exposures with organic reactions, however, these relationships have proven to be unclear or not reproducible. It has already been demonstrated in individual subjects that the variation in the rate of heart rate variability (HRV) decreases or the frequency band of the HRV under EMF exposure is narrower, with the reference value being a HRV control measurement before EMF exposure for comparison. Such a decrease in the range of HRV variation also occurs in the presence of any other stressor on the subjects and is the focus of many sports science and health prophylactic studies (high HRV is generally considered indicative of healthy physiological regulatory ability, whereas low HRV indicates cardiovascular health impairment System as well as the general organic constitution). Also, changes in the number of HRV fundamental harmonic waves were found under EMF influence, which, in combination with measurement of the capillary bed microcirculation by means of laser Doppler methods and the electrical potential distribution on the skin surface of the subject, indicate EMF sensitivity (eg Tüngler, A Von Klitzing, L, 2013: Hypothesis on how to measure electromagnetic hypersensitivity, Electromagnetic Biology and Medicine 32: 281-290.).
However, general or absolute or mean increases or decreases in brain or cardiac function over a test period do not provide a means of demonstrating a prior association between high-frequency EMF and such organic functional aberrations. So far no reliable criteria could be found, which indicate a similar causal relationship. The above methods for EEG or ECG / HRV-based analysis of vegetative regulatory dynamics under EMF influence are subject to many subjective influences on the part of the subject under test and therefore do not provide a reliable indication of the presence or absence of EMF sensitivity.
The object of the present invention is to individually determine the indication of an EMF exposure reduction by being able to detect human biologically relevant interactions of high-frequency EMF with the human organism in a rapid and meaningful manner. The association of specific radiofrequency EMF immission trajectories with autonomic regulatory disturbances should be evidently detected. Costly differential diagnostic measures and laboratory analyzes to determine an EMF sensitivity should be reduced. A corresponding test and measurement arrangement should be able to be carried out in particular also by untrained person locals independently and in an uncomplicated way.
Since in the present method there always takes place a consideration of a test subject or its location, ie immission-side viewing of high-frequency electromagnetic fields, these are referred to in the following as high-frequency immissions (HFI). However, since HFI and high frequency EMF always emanate from a technical emission source, the same could be said of high frequency emissions (which, however, decrease with increasing distance to the emission source).
The above objects are achieved by a method having the characterizing features of claims 1 and 2.
If the following is a description of a course of the high-frequency immissions (HFI) or of heart rate variability functional parameters (HRV-FP), then in each case the changing height or intensity, considered over a defined time interval, of one of the HFI, HRV-FP These values are measurement and / or calculation quantities that are based on measurement and calculation methods that are used by standardized measuring devices or processor devices for determining HFI, as well as other ECG-based measuring devices or processor devices for determining the HRV-FP are available according to the prior art (see for example DE 603 06 856 T2, DE 102006 039 957 B4).
The heart rate variability determination used in the course of the analysis method according to the invention is a recognized medical-technical method for evaluating the physiological state of a test subject. In the field of HRV, there already exists an extensive state of the art or a multiplicity of patent applications relating to specific applications and modified methods. DE 102009002134 A1 discloses, for example, a method for determining the effects of a magnetic field on a test subject, wherein different HRV (total) values can be determined before and after a change in the magnetic field has been made.
In general, the HFI value (in each case related to a specific instant) designates the power flux density or an equivalent measured variable for representing the intensity of electromagnetic high-frequency immissions acting on the test subject, which can be detected by means of a dosimeter or spectrum analyzer, which is further substantiated in the bottom, together with suitable reception subscribers are. Here, the unit [ower / area] is common, e.g. mW / m2 or μ W / m2. However, high frequency immunity levels can also be expressed in alternative units, which in most cases can be interconverted by simple mathematical operations. In rare cases, for example, data in V / m, A / m or in logarithmic (dB) units will be available.
The HRV-FP values, which are also substantiated in the following description, each denote performance values of vegetatively controlled organ functions of the test subject such as e.g. Cardiac and Respiratory Rhythms, Ratio of Sympathetic / Parasympathetic Activation, Respiratory Sinus Arthtenia, SDNNrr. Regardless of their organ or rule-specific statement about the organism of the test subject and their specific, related to specific body rhythms unit the above parameters are based on an analysis of the variance of the time interval of determined by Elektrokardio¬GG measuring method (ECG), successive heartbeats. In principle, all HRV FPs may be considered "Time Domain Measures" or "Frequency Domain Measures". known sowieadaptierte parameters are used, which are calculated from standard provided by an ECG measuring device RR interval sequence data.
The measurement of the heartbeats and ECG determination takes place in a manner known per se by means of pulse probes applied to the test subject. An objective analysis of the HRV is performed today nachstandardisierten mathematical methods, wherein measurement data of the heart rate (fundamental frequency) are transmitted, for example by means of Fourier transformation of the time domain in the frequency domain and can be displayed as a power spectrum. The representation of this Leistungsspekt¬rums is referred to as a spectrogram or AutoChrone image. Based on the data contained therein, according to the prior art various other frequencies present in the heart rate may be visualized that have been empirically assigned to other body rhythms (e.g., breath, blood pressure, blood flow, sympathetic, parasympathetic).
The HFI and HRV-FP traces are graph-representable in a two-dimensional coordinate system, but may also be described purely by mathematical functions, algorithms, and processor-usable data sets, i. be resolved alphanumerically. Regardless of their presentation or graphic visualization, the HFI and HRV-FP values are each a function of time, i. Each time point, which can preferably be represented on an abscissa of a diagram, is assigned a value that can preferably be represented on an ordinate of the diagram. Of course, the assignment of abscissa and ordinate can also be done in the opposite way, i. the time course also in the vertical direction or in the form of a third spatial axis of a three-dimensional diagram. The time frame for acquiring and / or displaying the HFI, HRV-FP values may be arbitrarily selected according to respective analysis requirements and processor capacities; it can e.g. from a few milliseconds [ms] to several seconds [s].
The two-dimensional diagram form for representing the indicated courses can also be supplemented by a third dimension in order to provide additional information about respective measurement or calculation quantities, which are shown in the graphical representation e.g. be expressed by a corresponding coloring, shading or hatching. In any case, a classic diagram with a horizontal X-axis as the time axis and the vertical Y-axis for mapping the HFI, HRV-FP values is preferred. The HFI, HRV-FP profiles resulting from the HFI, HRV-FP values resulting in a reference period are thus preferably subjected to a further evaluation as function graphs.
It should be noted that the indicated courses for carrying out the method according to the invention do not necessarily have to be represented in diagram form or in the form of a graph, but can also be analyzed and compared in a purely computational manner. Also, the result of an assessment of the presence of human biological relevant influence of the test subject could thus also be in purely alphanumeric form or by another, e.g. visualization represented by colors and / or symbols.
According to claim 1, the following method steps are provided: Determining the course of an intensity of radio-frequency immissions (HFI), preferably measured as power density density, with respect to the test subject over a defined period of time Determining the variance of the time interval of the electrocardiogram measuring method (ECG) , successive heartbeats and the course of determined therefrom, associated with the heart rate variability (HRV) vegetative function parameters (HRV-FP) of the test subject over a defined period (from the prior art, several, respective organ districts associated vegetative function parameters or history diagrams are known, which from an objective detection According to standardized evaluation software of the HRV technique, it is possible, from a frequency analysis of the temporal variance of successive heart beats, to have further orgasm present in the heart rate to determine niche frequencies. These frequencies calculated from the heart rate using the Fourier transform mathematical method were empirically assigned to specific body rhythms or vegetative functions such as sympathetic, parasympathetic, blood pressure, respiratory, etc. The determined HRV-FP traces may optionally be represented in the form of a graphical history, alphanumeric or algorithmic data, functions or function graphs and, if necessary, e.g. to graphically visualize the result of the analysis, to be transformed into an alternative representation.)
Frequency ranges or HRV-associated functional parameters determined from the HRV correspond in each case to a specific vegetative function such as, for example, Blood pressure, blood circulation, respiratory activity (depth and frequency), sympathetic, parasympathetic, thermoregulation, other organ activities; the assignment of HRV-extractable frequencies to respective autonomic functions of the human organism has been empirically determined and is detailed in papers on heart rate and HRV detection in accordance with the prior art. The HFI and HRV-FP profiles are synchronized at least in sections and calibrated for at least one common reference time interval (it is also possible to compare a plurality of reference time intervals or sub-time intervals which are congruent with respect to the chronology, the actual synchronization and calibration can either in hindsight or after data export of respectively inserted HRV and HFI analyzers or also in real time or by means of a data link of the HRV and HFI analyzers.) - performing an analysis of the HFI history in which primary significances, ie temporary, within the reference time interval, time-limited deviations of the HFI curve from defined static or dynamic mean values or HFI reference data, in particular significant increases and / or decreases in the HFI curve with respect to temporally preceding and / or subsequent sections / phases of the HFI curve are detected Any significant deviation from the usual frequency spectrum of HFI measurements is considered as HFI reference data or control curves, either average values from high frequency intensities or power flux densities determined at the site of the test subject or in environmental medical recommendations such as the precautionary evaluator Salzburg Chamber of Physicians or the "SBM standards of building biology metrology" defined high-frequency intensities and power flux densities.) - Performing an analysis of the HRV-FP course, in we secondary significance, i. temporary, i.e. Within the reference time interval, time-limited deviations of the HRV-FP curve from HRV-FP reference data or control courses, which are drawn from empirical data and / or generated dynamically, i. are derived from current measurement / calculation data or HRV-FP progressions or extrapolations, in particular significant increases and / or decreases in the HRV-FP progression compared to chronologically preceding and / or subsequent sections / phases of the HRV-FP progression. In this case, in turn, any deviation which exceeds the technical standard given oscillation width at HRV-FP measurement or calculation values is regarded as significant. [Note: in any case, the deviations or signatures to be detected within the reference time interval last less than the total reference time interval considered, i. are detectable as a limited temporal event or as a swelling and / or decongesting phenomenon. In contrast to known biofeedback methods, the event history of the respective HRV-FP profiles is thus analyzed. In the special case of a shielding experiment, described in more detail below, it might prove sufficient to detect only a swelling or waning of a secondary significance in a given time interval or at some point in time as a sufficient criterion for an analysis score; the same also applies to the detection capability of the primary significances or temporary deviations of the HFI curve]
A high detection accuracy of the analysis method according to the invention can in particular in the case of generation of HRV-FP reference data or runs using a combination of empirical data or in the past from measurement results on the same test subject or third test subjects received rule courses with based on the current test of the test subjects determined HRV -FP measurement / calculation data. HRV-FP and / or HFI reference data can either be drawn from empirical data or generated dynamically, i. derived from current measurement / calculation data or HRV-FP and / or HFI courses or extrapolated. The reference to HRV-FP and HFI control curves can be made by defining a defined fluctuation range or a defined frequency of irregularities of the courses, which are not yet to be evaluated as significances. By depositing such analysis algorithms Fehlinterpretati¬onen can be largely avoided and the analysis result can be optimized.
Upon detection of a defined number of coincident, i. essentially simultaneously occurring primary and secondary significances in the HFI course and in the HRV-FP course, there is a positive assessment of the presence of a human biologically relevant influencing of the vegetative regulatory mechanisms of the test subject.
Thus, there is a clear temporal assignment of secondary significances detected in the HRV-FP curve to primary significances detected in the HFI curve. The occurrence of secondary signatures in the temporal HRV-FP course is causally related to the occurrence of primary significance in the temporal HFI course.
As a consequence, exposure of the test subject to those HFI sources which have been determined to be relevant by the method of the invention can be reduced. In the simplest case, such an exposure reduction can be done by changing the sleeping place or by changing from wireless to corded home and communication technology. Exposure reduction is also possible by technical shielding measures on the building substance, e.g. using grounded plaster / wire mesh structures or carbon fiber coatings, such as those offered by the building materials producer Sto or Rigips especially for HFI reduction.
The detection of the coincident significances preferably takes place automatically, but it can also be done in a manual way e.g. based on a comparison of graphical visualizations temporally with one another corresponding HFI and HRV-FP courses. By means of a separate graphical representation of the correlation of HFI and HRV-FP courses, the relationships can be clearly displayed. The HFI and HRV-FP profiles can hereby be arranged one above the other in diagram form, with one another or one above the other. In a preferred variant of the method, the graphic output of the HFI and HRV-FP profiles already takes place in such a form of representation calibrated with respect to the analysis period or reference time interval. To quickly associate significances in the HFI history with corresponding significances in the HRV-FP history, the graphic representations of the two courses are arranged exactly one above the other, wherein the graphical representation of the HFI progression can have vertical grid lines, each of which characterizes a specific time point of the analysis period. These raster lines are aligned either with corresponding raster lines indicating the same point in time in the graphical representation of the HRV-FP curve or run continuously with the latter raster lines. The raster lines are at a suitable time interval of e.g. 10-30 minutes, preferably equidistant from each other. An optional representation would be if the HFI curve and HRV-FP curves are mapped in different colors and overlapping each other.
Alternatively or in combination for determining an immediate correlation of a trigger from the HFI with HRV-associated autonomic functional parameters (HRV-FP), the determination of such a correlation can also be carried out indirectly: Determining the variance of the time interval by means of electrocardiogram measurement methods ECG), successive heartbeats and the course of determined therefrom, associated with the heart rate variability (HRV) vegetative function parameters (HRV-FP) of the test subject over a defined period, - Comparison of the currently determined HRV-FP course with characteristic, HFl-induced secondary significance or HRV-FP reference data or rule progressions, which were respectively determined from the current test arrangement or from at least one of these preceding test arrangements with the current test subject or third test subjects and are now preferably in the form of digital or analog data, preferably in the form of graphic progressions, alphanumeric or algorithmic data, functions or function graphs on a storage device;
Note: The mentioned types of representations can optionally be converted into each other, e.g. for graphical visualization in graph form. The HRV-FP profile currently determined on the test subject can also be used in the form of a graphical progression, alphanumeric or algorithmic data, functions or function graphs, or converted into an alternative representation in each case. There is thus a definition of characteristic secondary significances occurring in the HRV-FP course, e.g. in the form of characteristic courses or deviations induced by primary significances in the HFI course of a corresponding experimental set-up. Secondary significances or characteristic courses or deviations of the HRV-FP course are preferably used here, which are recognized as a substantially coincident sequence of primary significances occurring in the HFI course during a corresponding test arrangement. However, it is also possible to use secondary significances or characteristic courses or deviations which occur in the HRV-FP course as a time-delayed sequence of primary significances recognized in the HFI course but are causally assignable. - Upon detection of a defined, at least partially present agreement of the current-determined HRV-FP course with characteristic, HFl-induced secondary significance or deviation of the currently determined HRV-FP course of the HRV-FP reference data, a positive assessment of the presence of a human biologically relevant influence the vegetative Regulati¬onsmechanismen the test subject takes place.
Note: In this case, each correlation between a significance detected in the event history of the HRV-FP profile and an HFI-typical significance stored on the memory device is detected as a match. The detection criterion of correspondence can also be defined negatively by defining HRV-FP reference data or control profiles and all significant deviations of the currently determined HRV-FP profile (both increases and decreases) from these HRV-FP reference data or control profiles Anomaly or secondary significance.
If appropriate, correlations which are defined in terms of HF1 may also occur with time-delayed, causative HFI loads or primary significances and be recognized as such.
The advantage of this embodiment variant according to the invention lies in the fact that in the course of the testing of a test subject, only one HRV-FP profile has to be determined from this and no HFI profile, thus no dosimeter, is required. The costs for a corresponding analysis can thereby be reduced, especially since measurement arrangements or analysis processes with regard to high-frequency immission (HFI) are only required for generating basic analysis algorithms with respect to characteristic correlations with HRV-FP profiles, but not with all the following, the customer or customer Analysis subsystems provided with the test object, in which only an HRV analysis device office of associated storage unit with analysis algorithms deposited thereon is to be provided. Any inaccuracies in the result, in particular the analytical method simplified according to the measuring device arrangement according to claim 2, can be compensated for by the following measures for increasing the accuracy of the analysis method.
It should further be noted that it is possible to modify the below-mentioned, associated with heart rate variability (HRV) vegetative function parameters (HRV-FP) by mathematical methods or specific frequency analysis method or transform into other measures and units, without departing from the inventive concept proposed analysis ,
A determination of characteristic significances, courses or deviations in the HRV-FP course can take place on the basis of third reference test subjects in a standardized analysis arrangement and then serves for comparison with HRV-FP courses of subsequently test subjects to be analyzed. The characteristic HRV-FP significances, courses or deviations can also be determined directly from the test subject to be analyzed, e.g. during a trial exposure to one or more HFI sources, described in more detail below.
The characteristic HRV-FP significances, courses or deviations can be stored on a storage device both as a static data record and as a dynamic data record which can be continually expanded or optimized on the basis of new analysis results. The characteristic HRV-FP significances, courses or deviations or their detection prevalence used for a comparison according to the invention can also be modified as a function of individual parameters of the current test subject. As individual parameters, factors related to both the physical constitution of the test subject (eg individual ability to regulate, medical diagnoses, examination results, certificates of the presence of an edrosensitivity, environmental medical reports, etc.) can be used as well as with the environment of the test subject or with the person living in it or work environment or exposure to radio frequency emissions detected or expected at the location of the current measurement arrangement (eg, type, number, distance of adjacent HFI sources, frequency ranges, timing, transmit powers, power flux densities, etc.) By reference to such, related to the test subject an improved analysis result can be achieved with individual parameters, the probability of overlooking a human biologically relevant influence by unrecognized HFI sources or due to individual compensation mechanisms of the test subject is reduced.
The characteristic significances, courses or deviations in the HRV-FP profile can be stored in the form of digital or analog data, in particular in the form of graphic progressions, alphanumeric or algorithmic data, functions or function graphs on a memory device in data communication with the comparator.
In other words, in both method variants presented above, an analysis of HRV-FP progressions is always carried out with regard to the presence of section-wise or temporary congruences with defined (HRV-FP and / or HFI) reference characteristics or courses. The congruence or secondary significance here does not necessarily have to be determined graphically, but can also be detected in a purely computational manner. Congruences to be detected can each be defined by defining permissible fluctuation tolerances and similarity criteria. Optionally, a congruence or secondary significance may also occur with time delay to a causative, HFI-related event and be detected as such.
In the present exemplary embodiments according to FIGS. 1-29, the most significant HRV-based autonomic functional parameters (HRV-FP) were used for the analysis according to the invention. It is possible to derive further possibly refined vegetative function parameters from the HRV basic evaluation or the HRV-FP, in which again the described significances are reflected, without deviating from the idea according to the invention.
According to a first preferred embodiment variant of the analysis method according to the invention, a heart rate variability (HRV) frequency information record is analyzed as the vegetative function parameter (HRV-FP) (in the case of graphic visualization, this is also referred to as "spectrogram" or "auto-chrone image") the over a defined Frequenz¬ bandwidth from substantially 0 to 0.5 Hz reaching activation each a specific
Frequency range of associated vegetative organ functions of the test subject represents / reproduces, in particular: 0.04 to 0.15 Hz: Low Frequency (LF) with correspondence: predominantly sympathetic activity, to a lesser degree also vagus activity, assignment in particular of the blood pressure and circulation rhythm (note: the vagus is the the largest nerve of the parasympathetic nervous system and involved in the regulation of the activity of almost all internal organs, the terms parasympathetic and vagus are used interchangeably below), 0.15 to 0.40 Hz: High Frequency (HF) with correspondence: vagus activity; Assignment of, in particular, respiratory functions, preferably the respiratory sinus arrhythmia (RSA) reflecting the modulation of the cardiac rhythm by respiration, and activating respective frequency ranges or autonomic organ functions from the HRV frequency information record in the form of amplitude intensities preferably visualized by color coding, wherein as HF1 induced secondary significance is detected when spontaneous activation of previously uneffected or negligible intensity regions or vegetative organ functions over a bandwidth (each understood as delta) of at least 0.05 Hz, preferably at least 0.1 Hz, more preferably across a band width of more than 0.2 Hz. Activation is understood in the present context in particular when the amplitude strength of respective frequency ranges or vegetative organ functions increases by more than 20%, preferably by more than 50%.
The HRV frequency information record or spectrogram is a clear representation of complex rhythm information contained in heart rate and heart rate variability, respectively. Here, the information i.d.R. represented in three dimensions: abscissa = time, ordinate = frequency, color = amplitude / strength of activation of respective vegetative functions; however, the subject information may also be evaluated as processor-usable alphanumeric code).
Frequencies represented within the subject frequency bandwidth arise from the modulation of the heartbeat, i. by changing the time intervals between directly successive heart beats (= RR intervals). In addition to the heartbeat rate calculated as the average of beats per minute, various other rhythms can be detected in the heartbeat used to control other organ systems. In the case of the coupling between cardiac and respiratory rhythm (QPA) during sleep, described below, the arrhythmia is transmitted to the heart rate and thus becomes visible in the spectrogram (RSA). The information content given here is of great importance for the assessment of the sleep architecture or its chaotisation by HFI loads.
In the method variant described above, it is detected according to a preferred Auswertever¬fahren as secondary significance, if a short-term, preferably not more than ten minutes, more preferably not lasting more than a minute activation of previously nichtbzw. in negligible intensity of activated frequency regions or vegetative organ functions, which in the case of a graphical visualization of the HRV frequency information data set represents an approximately needle-shaped elevation running orthogonal to the time axis. In a preferred variant of the method, a sufficient criterion for the detection as a secondary significance must be a simultaneous activation of more than 50%, preferably more than 70%, of the entire frequency band of the HRV frequency information data set or spectrogram comprising essentially 0 to 0.5 Hz.
In the evaluation of the HRV frequency information data set or spectrogram, it can also be evaluated as activation or as secondary significance if the amplitude strength of respective frequency ranges or vegetative organ functions increases by more than 20%, preferably by more than 50%, and the analysis method as well may include the possibility of weighting detected secondary significances on the basis of respectively determined color-coded amplitude strengths in the case of graphical visualization and / or on the size of the respective activated frequency bandwidth, wherein secondary significances at which high amplitude intensities are detected. have been detected over a wide range of the frequency bandwidth of the HRV frequency information data set activating activations are subjected to a larger factor for the assessment of the presence of human biologically relevant influence of the vegetative test subject regulation mechanisms than those secondary significances, at which lower amplitude intensities or over small ranges of the frequency bandwidth of the HRV frequency information record were detected activations. By evaluating events or secondary significances according to their intensity and thus their potential impact on the autonomic regulatory system of the body, the expressiveness of the method according to the invention can be substantially increased. The performance of such weighting as a function of the intensity of the events / secondary significances can also be carried out analogously in all other HRV-FP processes mentioned.
In a special analysis method, in turn, the respiratory sinus arrhythmia (RSA), which can be represented in the HRV frequency information data record or spectrogram, or its course lying in the range from 0.2 to 0.3 Hz, is brought into focus, with an upper end of the frequency bandwidth substantially convexly curved course of activated organ functions is assessed as HFI-induced secondary significance, the convexity preferably extending over a time interval of between 30 and 160 minutes, more preferably between 60 and 120 minutes.
The already mentioned respiratory sinus arrhythmia (RSA) is in each case the high-frequency variability of the heart rate, which reflects the intensity of the modulation of the heart rhythm by the respiration and can be determined by standardized HRV calculation methods. The RSA is converted to the logRSA by means of a decade logarithm, which according to a further variant of the method of analysis according to the invention can also be investigated as autonomous autonomous functional parameter (HRV-FP) with respect to HFI-induced significances.
Note: The course of each organ function analyzed in the HRV frequency information data set or in the spectrogram is not necessarily shown in the spectrogram as a solid line, as with other vegetative HRV function parameter curves described above, but as a more or less fragmented area cloud, with each Point or each area of this cloud is associated with information about the amplitude strengths of respective organ systems at certain times. An objectively described convexity of the RSA profile with respect to the time axis is generally preceded by an attenuation of the amplitude intensity under the influence of HFI, which in the case of a graphic visualization or in the spectrogram, for example, as a paler or more directional blue spectrum from previously even more red / yellow spectrum-tending coloring. Such an attenuation of organ functions corresponding to the respiratory sinus arrhythmia (RSA) in the frequency range of substantially 0.2 to 0.3 Hz of the HRV frequency information data record can additionally or alternatively already be evaluated as HFI-induced secondary significance in the method according to the invention.
It has surprisingly been found in the course of the development of the method according to the invention that the above-mentioned time interval (30 to 160 min or 60 to 120 min), in which the organ functions lying in the frequency range of about 0.2 to 0.3 Hz with respect to their frequency abschwellen or show a convex course, with the so-called BRAC (basic resting
Activity cycle) corresponds. The BRAC designates an ultradian chronobiological rhythm with a period of approximately two hours, comprising an activation phase of i.d.R. 80 to 120 minutes and a regeneration phase of i.d.R. 20 to 30 minutes. In sleep, this rhythm manifests itself as a change between deep sleep and REM phases, the latter manifesting themselves at the lowered arc shapes or end portions of the convex RSA progress sections. While during the day a rhythmic change between about 90 min activation or sympathetic activity and about 30 min regeneration or parasympathetic / vagus activity shows (Note: the periodic90 / 30 min change of BRAC is an average time, which this ratio may be reversed at night or during sleep: the parasympathetic / vagus activity is about 90 minutes, the sympathetic activity about 30 minutes.
The BRAC-characteristic course of RSA activation in the HRV frequency information record reflects activation by the sympathetic and parasympathetic (rather than the respiratory center). A reflection of the BRAC in the otherwise substantially bar-shaped or parallel to the time axis (in a frequency range of approximately 0.25 Hz) extending RSA course during sleep is to be regarded as a pathological tendency and thus represents a suitable criterion to a HFl- to detect induced disturbance of the vegetative regulatory mechanisms of the test subject.
An increase in RSA progression associated with weakening of the respiratory center, e.g. at approximately 0.3 Hz indicates an anti-cyclical slowing of respiration
Decreasing the RSA course indicates an anticyclic acceleration of the test subject's breath and thus increasing sympathetic activation; in total, the respiratory rhythm is at least partially chaotized.
According to a further embodiment variant of the analysis method according to the invention, the temporal course of the heart rate of the test subject is analyzed as the vegetative function parameter (HRV-FP), and the occurrence of short-term arrhythmias in relation to the chronologically preceding and following heart rate course is detected as a secondary significance.
In a preferred variant of the analysis, the arrhythmia of the heart rate and thus secondary significance (1'-n ') is the occurrence of extrasystoles, whereby an extrasystole over the course of the heart rate measured in the unit: [number of heartbeats per unit time, preferably per minute] is considered to be short-term , represents a sudden increase in the heart rate over the respectively preceding and subsequent heart rate progression and thus as a significant shortening of the interval duration between two successive heartbeats. In a preferred variant of the method, sudden increases in the heart rate are recognized as extrasystoles if these increases are in each case more than 30%, preferably more than 50%, particularly preferably more than 100% of the immediately preceding or following heart rate or the average heart rate determined in the respective measurement period.
Colloquially, an extrasystole is also referred to as a double stroke or intermediate heart beat. In the analysis method according to the invention, in particular intermediate beats are detected in which the heartbeat-inducing stimulus does not originate from the sinus node of the electrical cardiac conduction system, but is to be considered as a consequence of high-frequency immissions. in the case of a graphic presentation or evaluation, extrasystoles would map as needle-shaped elevations in the heart rate course. That portion of the heart rate progression corresponding to an extrasystole corresponds to a curve with a very strong slope approaching infinity. In practice, extrasystoles are considered to be approximately orthogonal to the time axis (thus, i.c., vertically to the abscissa of an HRV plot). The greater the short-term rash of heart rate. a corresponding curve, the shorter the time interval between two consecutive heartbeats of the test subject.
According to a further preferred embodiment variant of the analysis method according to the invention, the temporal course of the respiratory sinus arrhythmia (RSA) reflecting the respiratory modulation of the heart rhythm is analyzed as a vegetative functional parameter (HRV-FP), which is preferably converted to logRSA by means of a decadal logarithm. As previously stated, respiratory sinus arrhythmia RSA is the high frequency variability of heart rate which reflects the rate of cardiac rhythm modulation by respiration. RSA or logRSA is also a measure of the tone of vagal activity.
Alternatively or additionally, the temporal course of the pulse-breath quotient (QPA) of the test subject can be used as the vegetative function parameter (HRV-FP), wherein preferably those events are detected as secondary significances in which the QPA course of a ratio of pulse / breath = 4: 1 deviates by more than 20% and / or in which temporary peaks or maximum values of the QPA gradient occur.
The vegetative function parameter of the QPA (pulse / breath quotient), which is also calculated in the standardized manner from the modulation of the heart rate, indicates how often the heart beats (= pulse waves) during a respiratory cycle. The QPA increases with tension and decreases with relaxation. The interplay of heartbeat and breath is referred to in chronobiology as the center of the rhythmic system in man. The organism seeks to keep the equilibrium balance between pulse and breath constantly in a rhythmic balance. During undisturbed sleep, optimal synchronization between pulse and breath is in a 4: 1 ratio. It has been found that this vegetative function parameter is particularly suitable for detecting an HF1-induced load of the test subject.
In a preferred embodiment variant of the analysis method, all peaks or maximum values occurring within the QPA profile are detected as secondary significances. These are in the case of a graphical representation of the QPA course again as approximately needle-shaped or serrated elevations. It is hereby preferably checked according to claim 1 whether coincident primary significances in the HFI profile are also present at those points in time of the reference time interval at which such QPA peaks were determined.
According to a further preferred embodiment variant of the analysis method according to the invention, the temporal course of the vegetative quotient (VQ) of the test subject is analyzed as the vegetative function parameter (HRV-FP), which is determined by the ratio of low frequencies assigned predominantly to the sympathetic, in particular to the blood pressure rhythm, substantially between 0.04. 0.15 Hz extending HRV frequency ranges (LF) to the parasympathetic nervous system, in particular the high-frequency frequency bands associated with the attenuation substantially between 0.15 - 0.40 Hz (HF), whereby secondary significances are detected: - temporary increases of the vegetative quotient (VQ) more than 20% (compared to immediately preceding VQ values or to a VQ mean value determined in the reference time interval) or / and excesses or sympathetic dominance of the vegetative quotient of more than 3: 1, preferably more than 5: 1 ,
The HF (high frequency) range comprises vegetative rhythm fluctuations with period lengths of approx. 2.5 seconds to 7 seconds, or 0.15 - 0.40 Hz, respectively. HF power performance corresponds to the activity of the parasympathetic nervous system and mainly reflects heart rate variations due to modulation via the respiration. The LF (low frequency) range covers the frequency range of about 7 -25 seconds or 0.04 - 0.15 Hz). Performance in this frequency range is primarily affected by the sympathetic nervous system and, to a lesser extent, the parasympathetic nervous system. This frequency range was formerly also called baroreceptor area, since the activity of this receptor is very reflected here. The low frequency components of heart rate variability correspond in particular to the blood pressure rhythm.
It is again preferably checked according to claim 1 whether coincident primary significances in the HFI sequence are also present at those times of the reference time interval at which such elevations or sympathetic dominances were ascertained.
According to a further preferred embodiment variant of the analysis method according to the invention, the time profile of the SDNNRR (standard deviation of normal-to-normal intervals) of the test subject is analyzed as vegetative function parameter (HRV-FP), i.e. an HRV-associated statistical spread measure about the mean of the heart rate interval duration or its differences representing the change in the total variability of substantially all frequency ranges of artifact-adjusted RR heartbeat interval series within a specified time interval. The DNNrr is a measure of overall variability over all frequency ranges. Evaluation as HFI-induced significance occurs especially in segments with spontaneously elevated SDNNRR, but can also be evaluated according to other criteria described below. Instead of the SDNNRR, a corresponding or modified statistical characteristic ("time domain measure") such as e.g. SDNNIDX / ASDNN of the test subject.
The activation of autonomous regulatory mechanisms or vegetative function parameters during the daytime or in the wake of a multitude of influences, which are directly related to physical or mental activities of the test subject. Therefore, by determining the HRV-associated vegetative function parameters in a time interval and using them for an analysis according to the invention with regard to HF1-induced significances in which the test subject sleeps, influences attributable to the subjective sphere of the test subject can be largely eliminated. Since extensive activation of parasympathetic entry and external variables are restricted or reduced during sleep, immediate response of the autonomic system to HFI immissions can be detected with significantly greater evidence than is possible with two-day biofeedback procedures. The requirement for testing during sleep is particularly indicated in the case of an analysis of respiratory sinus arrhythmia (RSA), the pulse-to-breath ratio (QPA) or the SDNNRR-related vegetative function parameters (HRV-FP), in particular This underlying synchronization of the heart rhythm by breathing usually only occurs during sleep. Sleep during sleep is also indicated in an analysis of the activation of other organ systems in the HRV-frequency information dataset or spectrogram.
In a particularly preferred variant of the method, all events or course segments are detected as secondary significances in one of the above-mentioned HRV-FP courses, which graphically shows a gradient k> 1 or a pitch angle α of more than 45 °, preferably of more than 70 ° , Particularly preferably of approximately 90 ° relative to the time axis (slope k = °°) or an approximately needle-shaped elevation of the HRV-FP curve, wherein both aforementioned events or detection criteria each exceeding the average HRV-FP curve must be at least 20%, preferably at least 30%. As an average HRV-FP curve or as a detection reference value, either a HRV-FP value averaged over the entire period of the test or the reference time interval can be used, or a HRV-FP value averaged over a shorter period, e.g. a detected event or secondary significance before and / or after a maximum period of 120, 60, 40, 20, 10 or 5 minutes. The aforesaid slope or slope angle respectively applies to a curve tangent applied to the history graph at a particular time.
According to a further preferred variant of the method, it is provided that at least one sub-time interval is opened from the time of detection of HFl-induced secondary significances, wherein at least one further, within the range of the positive evaluation of the presence of a human biological relevant influencing of the vegetative regulatory mechanisms of the test subject Sub-time interval occurring secondary significance must be detected and whereby the opened sub-time interval is preferably a maximum of 5 minutes, more preferably a maximum of 1 hour.
In particular, sub-time intervals may also be determined according to the manner of the course of the power density of the HFI in the reference time interval, i. as a consequence of the detection of a primary significance in the HFI course.
A sufficient condition for the positive assessment of the presence of a relevant HFI-indicated influencing may furthermore be the requirement of detecting a defined number of secondary significances, preferably at least three, more preferably at least five, per reference time interval, wherein the reference time interval is preferably between 0.5 and 12 hours, more preferably between 5 and 30 minutes.
In practice, the measured values or data recorded by the ECG measuring device or HRV-FP evaluation device and the HFI measuring device are provided in different recording frequencies. In order to make computational comparability of the ECG or HRV-FP data and the HFI data possible, it is intended to synchronize these. To synchronize the HFI and the HRV-FP curve for the considered reference time interval (Τ ·, .eta.), An equal number of the HFI and HRV-FP measurement or calculation values constituting the two events are respectively generated, preferably by interpolation and / or Extrapolation and / or selection of acquired HFI and HRV-FP measurement / calculation values. Time-corresponding HFI and HRV-FP measurement / calculation values are each provided with the same time indices. In a preferred embodiment, the reference time interval for both (HFI and HRV-FP) courses is subdivided into a multiplicity of sub-time intervals of a maximum of 30 seconds, particularly preferably of a maximum of 10 seconds. Here, in each sub-time interval in each case that HFI measurement / calculation value and that HRV-FP measurement / calculation value is selected, which is the most significant with respect to the respectively defined as detection criterion for primary and secondary significances characteristic curve, preferably e.g. Minimum, maximum, median, average or delta values. The data obtained from two different signal channels (the HFI measuring device and the ECG measuring device or the HRV-FP evaluation device) are thus brought to the same number of measurement / calculation values by means of block operations.
As a sufficient condition for the positive assessment of the presence of a relevant HFl-indicated influencing, it can also be provided according to a further preferred variant of the method that at least one, preferably at least three coincidental significances (in the HFI and HRV-FP course) must be detected per hour or in the analyzed reference time interval, more than 20%, preferably more than 30%, particularly preferably more than 50% of the primary significances detected in the HFI curve correlate with coincidental secondary significances in the HRV-FP curve.
According to a preferred variant of the method, it is provided that one or more of the following characteristics are preferably detected as primary significances in the HFI curve: maximum values or power flow peaks, which in the case of a graphic representation can be seen as approximately needle-shaped or jagged amplitude excursions and / or - approximately staircase-shaped course sections, in this case preferably the beginning and / or end sections of a respectively approximately horizontal or a consistent power density density indexing course section of the staircase form and / or ebendiesen approximately horizontal run section itself, and / or - several, preferably more than three, preferably within a time interval of a maximum of 10, 30 or 60 minutes consecutive increases and decreases in the course (note: thus a comb-shaped waveform) and / or - exceeding a power flux density of more than 0 , 1 mW / m 2, preferably more than 0.05 mW / m 2, more preferably more than 0.01 mW / m 2 (Note: It is understood that the preceding mW / m 2 unit is also divided into other equivalent units for the determination of radio frequency Immissions or power flux densities can be converted) and / or - approximately constant power flux densities or in the case of a graphical representation substantially parallel to the time axis extending sections of the HFI curve and / or - reductions and / or increases in the HFI curve or the power flux density, which preferably have a difference of more than 10%, particularly preferably more than 30%, in relation to immediately preceding and / or following HFI progress sections, wherein as "immediately preceding" a period of several, eg 1-30 minutes, preferably from 0-60 seconds, more preferably from 0-10 seconds can be defined. (Note: The aforementioned detection criteria for increasing and / or decreasing HFI segments may also be used analogously to detect secondary significances in one of the HRV-FP methods.) Further, all of the aforementioned detection criteria may also be used to generate the characteristic, HFI-induced secondary significances HRV-FP reference data according to claim 2 are used.)
According to a specific variant of the method, it is provided that the interval duration, i. the time intervals of successive primary significances in the HFI course, and preferably also their intensity, i. the power density density difference between immediately preceding and / or subsequent HFI history sections is subjected to a frequency analysis and statistical values are generated for which respective thresholds or maxima and / or minima are determined, wherein exceeding or undershooting of these threshold values is recognized as the primary meta-significance and examining the occurrence of such a primary meta-significance on a temporal correlation with secondary significances detected in the HRV-FP course, the presence of such a correlation being considered an indication of the presence of a human biological relevant influence on the vegetative regulatory mechanisms of the test subject, and preferably as a statistical characteristic the SDNN (standard deviation of normal-to-normal intervals) or an equivalent time-domain characteristic related to the interval duration of successive primary significances will be.
In the course of the analysis method according to the invention, therefore, the course of the power flux density of the high-frequency immissions is not necessarily analyzed, but the HFI curve can be converted into statistical characteristics analogously to the HRV frequency analysis of the RR intervals discussed above, using mathematical functions such as Fourier transformation, which in turn are used as comparison variables or serve as primary significances for comparison with the occurrence of secondary significances in the HRV-FP course. For example, can be directly determined from the magnitude of these statistical characteristics such as (HFL-related) SDNN as well as derived and associated characteristics whether or not an HFI influence measured in a particular time interval is considered to be relevant or particularly burdensome. The classification of the relevance results here from the specific progression characteristic of the HFI curve measured by means of dosimeters, in particular from the number of LFD increases and / or LFD reductions and / or LFD peaks detected in a certain period and / or the ratio of each to successive ones Time-measured LFD intensities. The immission scenarios or primary significances described in the following exemplary embodiments can thus also be defined and detected in an indirect manner. For example, For example, a low or zero SDNN may indicate a less relevant HFI (primary) or meta-significance, while a relatively high SDNN indicates a particularly relevant HFI (primary) or meta-significance. The calculated statistical values or SDNN values are each provided with time indices for chronological assignability. Again, the aforementioned detection criteria can also be used to generate the characteristic, HFI-induced secondary significances or HRV-FP reference data according to claim 2.
According to a particularly meaningful variant of the method, it is provided that a targeted initiation of a high-frequency immission acting on the test subject is carried out by one or more dedicated test radio-frequency source (s) and thus one or more temporally determined primary significances are generated, wherein HRV currently determined -FP history of the test subject is checked for secondary significances coincident to the primary significance (s), and wherein the initiation of the high-frequency immission is preferably in a (related to high-frequency immissions, ie not caused by the test high-frequency source) HFI gradient range with approximately constant emission height over a defined period of time or in the presence of external high-frequency immissions which are below a defined threshold value, preferably below a power flux density of 0.4 mW / m 2.
According to a special variant of the method, it is provided that the test subject is protected during one or more time intervals by a shielding device preferably surrounding the test subject on all sides, such as e.g. is shielded from external high frequency immissions of the environment by electrically conductive metal mesh and the shielding device is either opened for a defined period of time and thus the test subject is exposed again to the external radio frequency emissions or the enclosed shielding device at the time interval of the shielding at least one test radio frequency source located within the shielding device is activated one or more times for a defined period of time, wherein the simultaneously determined HRV-FP curve of the test subject coincides with secondary significances occurring at the primary significance (s) generated by the test radio frequency source, and secondary, respectively, during the time of opening of the shielding device Significances is checked. If noticeable differences occur in the HRV-FP course of the test subject during the objectively artificial immission conditions (duration of the shielding and / or duration of the activation of the test radio-frequency source) compared to previous or subsequent time intervals, each with a comparatively regular immission level. High frequency power flux density), this is assessed as having a human biologically relevant effect on the vegetative regulatory mechanisms of the test subject.
By targeted debugging arrangements according to the preceding claims, the Untersu¬chungsdauer or the measurement operating time for the HFI and / or HRV-FP detection devices can be significantly reduced and the validity of the analysis result can be increased. l.d.R. For example, an evident analysis result on the presence or absence of an HFI-induced regulatory disorder can be given to the test subject within an examination period of 30 minutes.
In a preferred embodiment of the invention, it is provided to carry out a frequency-specific analysis of the HFI and in this case to compare several HFI courses with respective HRV-FP courses. In procedural terms, such comparisons may take place either simultaneously or sequentially. Viewed time windows can be arbitrarily fragmented. Also, a preferably automated analysis of selected HFI frequency bands or their HFI progressions, including comparison with time-corresponding HRV-FP progressions, can only take place at those times or intervals at which one or more defined significances are detected during the respective HFI progressions. The subject comparisons may e.g. typical HFI frequency bands include GSM, UMTS, LTE, WLAN, DECT, Bluetooth and other transmission standards. assigned.
From the high-frequency immissions acting on the test subject, it is thus possible to identify a plurality of HFI sources, each having different frequency band ranges, and to determine corresponding HFI progressions (intensity versus time), each time related to the test subject. In a preferred embodiment, in the method according to the invention, the frequency bands can be adapted to allergy-transmitting standards such as GSM, UMTS, LTE, WLAN, DECT, Bluetooth and the like. each detected as isolated HFI waveforms and compared with time-corresponding HRV-FP curves. It is also possible to use only those frequency bands or HFI profiles for the representational comparisons for which a frequency-specific exceeding of defined high-frequency immission reference values has been determined on-site by measurement; a definition of such immission reference values may in turn be e.g. based on environmental medical recommendations, such as the precautionary values of the Salzburg Medical Association or on the "SBM standards of building biology measurement technology".
Furthermore, it is also possible to combine or separate selected frequency bands by summing or subtraction or by other mathematical functions in a common HFI history, i. filter out certain frequency bands and thus exclude them from the analysis.
Such a multi-layered, frequency-selective analysis method of HFI for the detection of significant correspondences with the HRV-FP curves is advantageous in that it has been found in environmental medical practice that electrosensitive persons react in a very individual and difficult to assess manner to each existing high frequency immissions. Some individuals do not show a response indication to emissions from certain frequency bands, while emissions from other frequency bands (or specific clocks / pulse rates) are given a high response indication and vice versa. Since humans today are exposed to a growing variety of high-frequency emitters, both externally and internally, it seems particularly important in terms of emission reduction to determine those transmission standards or HFI sources which lead to the most significant individual burdens ,
A frequency-selective analysis method of the HFI may also be performed to generate characteristic, HFI-induced secondary significances in a trial arrangement according to claim 2 prior to the current experimental setup.
According to a special evaluation variant, it is provided that measurement or calculation values constituting the HFI profile and the HRV-FP sequence are respectively converted to a (chronological) gradient sequence (Ghrv-fp and GHfi), which consists of a multiplicity of temporally successive ones Gradients, ie Assemble difference formations between in each case two temporally successive HFI or HRV-FP measurement or calculation values. The gradient sequences (GHRv-Fp and GHfi) are formed in the case of a graphical representation as a bar graph and are analyzed in accordance with the invention in terms of their correlation.
According to a preferred variant of the method, corresponding gradient gradients (Ghrv-fp and GHfi) are created, wherein the gradient sequence corresponding to the HRV-FP course (GHrV-fp) and the gradient sequence (GHfi) corresponding to the HFI course overlap one another or side by side, in this case preferably mirror-symmetrical to a parallel or congruent with the time axis extending mirror axis are arranged. Preferably, at least one of the gradient sequences (GHrv-fp and GHFi) is scaled in such a way that orthogonal to the time axis average elevations of the gradients are approximately equal in both graphical processes.
According to a further variant of the method, envelope curves are applied to the graphic sequences of the gradient sequences (GHrV-fp, GHfi) which are to be compared with one another. In this way, any correlation in the time course of the gradient sequences (GHrv-fp> GHfi) or the simultaneous occurrence of primary and secondary significances can be shown particularly clearly. Even significances with comparatively small fluctuations or short-term duration can be visualized in this way.
According to another variant of the method, a correlation value, preferably in percent, is determined between the gradient sequence (GHrv-fp) corresponding to the HRV-FP curve and the gradient sequence (GHfi) corresponding to the HFI curve (thus for the HFI and HRV-FP determined within the considered reference time interval gradients). If this correlation value has a respectively fixed reference correlation value of e.g. Exceeds 30%, or 0.3 or 40%, or 0.4, there is a positive assessment of the presence of a human biological relevant effect on the test subject.
It should be noted that a previously described gradient formation with arithmetic and / or graphical evaluation can also be carried out analogously for an analysis of HRV-FP profiles according to claim 2 and dependent method variants, i. also the HRV-FP reference data or control curves or the characteristic, HFl-induced secondary significances can be compared in the form of gradient sequences with the respectively currently determined HRV-FP profiles of the test subject.
In order to exclude imponderables which result from a consideration of individual HRV-FP, in preferred embodiment variants of the analysis method, an analysis of several HRV-FP profiles with respective detection of secondary significances takes place. Thus, it may be provided that as a verification of the analysis, a comparison of the analyzed HRV-FP curve with at least one, preferably with at least two of the HRV-FP curves mentioned in the preceding claims with regard to the presence of substantially the same time or causally correlating secondary significances is performed. Preferably, the analyzed HRV-FP profile and / or the additional HRV-FP profiles used are in each case compared with in a time-corresponding HFI curve or with primary significances detected therein according to claim 1.
According to a further variant of the method, it is provided that, as verification of the analysis, a comparison of the analyzed HRV-FP profile in at least two, preferably in at least three of the HRV-FP profiles or analysis types cited in the preceding claims, respectively, a defined number of secondary significances are detected got to.
According to a further variant of the method, it may be provided that, depending on the determined frequency and / or standard deviation (size) of the HFI-induced secondary significances, a preferably in percent, gradual evaluation of the probability of the presence of an incompatibility of high-frequency immissions. In this case, the standard deviation is understood to be the deviation from respective metrologically determined mean values or expected values, which take into account typical measured value maxima and minima. In the simplest sense, the standard deviation may be defined as the variance of respective measurement or calculation values by their expected value. Their basic determination from the square root of their variance can be supplemented by extending functional parameters. For example, in determining the standard deviation, process-specific factors and characteristics such as hysteresis, resonance, addition or multiplication effects, typical swelling and decay effects, and the like (both related to HRV waveforms and HFI waveforms) may be included or filtered out.
Alternatively or additionally, a gradual, preferably in a percentage, evaluation of the vegetative regulatory capability of the test subject relative to HFl-induced secondary significances can be performed on the basis of a comparison with HRV-FP comparative data, these comparison data being based on current or past measurements on the test subject itself or on experience-external HRV -FP test series are based.
A graded, preferably percentile, assessment of the subject's vegetative regulatory ability against HFI-induced secondary significance may also be based on the rate of HRV-FP detected by the organism of the test subject after detection of HFI-induced secondary significance in HRV-FP again to a standardized or individually calculated HRV-FP control curve, that returns to defined setpoint values, take place (DieSoll values can in turn be stored statically on a database or dynamically generated in real time).
Claims 29 and 30 are directed to arrangements for carrying out the aforementioned Analysver¬fahren.
According to claim 29, this arrangement comprises a measuring device for measuring the intensity of high-frequency immissions (HFI), which is preferably determined as the power flux density, and an electrocardiogram (EKG) measuring device including measuring electrodes for determining the temporal variation of successive heart rate intervals and the course thereof determined with the heart rate variability (HRV) associated vegetative function parameters (HRV-FP) of the test subject and a processor-controlled Auswerteinrichtung together with this data-connected storage device on which one or more evaluation algorithms are deposited for performing an analysis method according to one of claims 3 to 28, wherein preferably the ECG measuring device and / or the evaluating device and / or the HFI measuring device are combined in a device that can be applied to the test subject, or the evaluating device is positioned at an external location and by means of the ECG-bz w. HRV-FP and HFI measuring devices determined data sets to this external, voresweisweiserverbased evaluation are transmittable for carrying out the analysis method according to the invention.
According to claim 30, the arrangement comprises an electrocardiogram (ECG) measuring device including measuring electrodes for determining the temporal variation of successive heart rate intervals and the course of determined, associated with the heart rate variability (HRV) vegetative Funktionspa¬rameter (HRV-FP) of the test subject and a processor-controlled evaluation together with this data link-connected memory device on which characteristic, HFl-induced secondary significances or progress characteristics in the form of digital or analog data are stored and one or more evaluation algorithms for performing an analysis method according to one of claims 3 to 28, wherein preferably the ECG measuring device and / or the evaluation ¬richtung are combined in a device applicable to the test subject device or the evaluation is positioned at an external location and by means of ECG or HRV-FP and HFI Messinrichtungener mediated data sets to this external, preferably server-based evaluation device for carrying out the analysis method according to the invention or for comparison with the characteristic, HFl-induced secondary significances can be transmitted.
The external evaluation device is preferably designed as a server-based device, to which the data records acquired by the ECG and / or HFI measuring devices, e.g. be transmitted in a wireless manner or by means of an internet-based upload. The server-based evaluation of the data sets in accordance with the analysis methods proposed according to the invention can be carried out fully or semi-automatically or with or without the supervision of a person skilled in the art. A client or the test subject is thus given the opportunity to obtain information on the presence or absence of an HFI-induced stress of his vegetative organism within a very short time. The calculation of the HRV-FP historical data made on the basis of the ECG measurement data can take place both at the designated server and directly in the device that can be applied to the test subject. Likewise, the analysis method according to the invention could also be carried out directly in the device applicable to the test subject or in conjunction with a stationary PC provided with a data transmission interface and a corresponding analysis software.
The invention will now be explained in more detail with reference to exemplary embodiments. Note: The following figures or diagrams are each color in the original, were resolved in pure black and white puncturing graphics to meet standard requirements. An additional information content of the graphics given by original color codings, for example amplitude magnitudes in the spectrogram, is explained in the technical description.
Show it:
2 shows a representation of the positioning of ECG measuring electrodes on the test subject FIG. 3 shows a heart rate diagram 4 determined from the detected RR intervals, a frequency-specific power spectral density 5 of respective frequency ranges determined in the sequence FIG. 4 shows a HRV frequency information data set which is represented as a three-dimensional spectrogram and contains the spectral information according to the preceding FIGS. 1-5 FIG. 7 shows a high-frequency immission (HFI) course measured over a period of 8 h (FIG. HFI peaks marked as primary significances)
8 shows a course of the heart rate of a test subject with detected extrasystoles, which is determined over the same period of time as the HFI curve according to FIG. 7, shows a radiofrequency immission (HFI) course measured by means of a dosimeter over a period of 8 hours
10 shows a course of the HFI measurement according to FIG. 9, represented as a spectrogram HRV frequency information record of a test subject FIG. 11, a graph of the vegetative quotient VQ of the test subject (sympathetic / parasympathetic activity. FIG )
FIG. 12 shows an HFI curve measured by means of a dosimeter over a period of 8 h (HFI peaks marked as primary significances). FIG.
FIG. 13 shows a course of the pulse-respiratory quotient (QPA) of a test subject determined over the same period of time as the HFI curve according to FIG. 12, with a focus on increasing the QPA according to HFI-induced secondary significances
FIG. 14 shows the HFI course measured over the same period of time as in FIG. 12 (HFI stair patterns or panels marked as primary significances).
FIG. 15 shows a course of the pulse-breath quotient (QPA) of the test subject determined over the same period of time as the HFI curve according to FIG. 14, with a focus on secondary significances of the QPAFig.16 corresponding to the staircase patterns of the HFI curve, measured by means of a dosimeter HFI course over a period of 8 h with primary signatures
FIG. 17 shows a course of the heart rate of a test subject with secondary significances, determined over the period according to FIGS. 16 and 17 over the same period of time as the HFI curve according to FIG. 16, of a test subject with a secondary test Significance Fig.19 shows an HFI curve measured by dosimeter over a period of 8 h with primary signatures
FIG. 20 shows a course of the heart rate of a test subject with secondary significances FIG. 21, determined over the period in accordance with FIGS. 19 and 20 over the same period as the HFI profile according to FIG. 19, of the SDNNrr (standard deviation of normal-to-normal intervals ) of a test subject with secondary significances
FIG. 22 shows an HRV frequency information record represented as a spectrogram during a sleep phase of the test subject
FIG. 23 shows a progression of the vegetative quotient VQ of the test subject (sympathetic / parasympathetic activity) shown over the period according to FIG.
FIG. 24 shows a progression of the heart rate of the test subject with the REM phases drawn in over the period according to FIG
FIG. 25 shows an HRV frequency information record represented as a spectrogram during a further sleep phase of the test subject (two nights after the patient) shown in FIG Analysis according to Figs. 22-25)
FIG. 27 shows a course of the vegetative quotient VQ of the test subject (sympathetic / parasympathetic activity) shown over the period according to FIG.
FIG. 28 shows a progression of the heart rate of the test subject with marked REM phases over the time period according to FIG. 26
29 shows a curve of the pulse-breathing quotient (QPA) of the test subject shown with reference to the preceding FIG. 26 with the REM phase Fig.30 drawn in, and an HFI curve measured by means of a dosimeter
FIG. 32 shows gradient sequences corresponding to the HFI and heart rate curves according to FIGS. 30 and 32, mirrored about a horizontal axis
FIG. 33 shows a heart rate progression 34 determined over the same time period as the HFI curve according to FIG. 30, the gradient sequences from FIG. 31 with applied envelope / trend curves, mirrored about a horizontal axis
FIG. 34 shows the gradient sequences from FIG. 33 with applied envelope / trend curves superimposed on one another
Fig. 1 shows RR intervals of successive heartbeats of a test subject taken in an electrocardiogram (ECG). A RR interval here refers to the distance between two R waves of an electrocardiogram measured along the time axis (here: abscissa) and thus the time interval between two successive pulse waves or contractions of the heart muscle of a test subject. For a frequency analysis described below, instead of the R-Zackenbzw. RR intervals other Interbeat intervals or derived indication sizes are used.
For exact measurement of the RR intervals, three ECG measuring electrodes E1 (negative electrode), E2 (positive electrode) and E3 (neutral electrode) are applied in a configuration shown in FIG. 2 at the breast area of the test subject. The detected sequence of RR intervals is referred to as a time series. In the present exemplary embodiment, the RR intervals were converted from a binary format into ASCII values and read in by an evaluation device. From the temporal variability of the RR intervals, the heart rate variability (HRV) can be determined by means of standardized mathematical operations. Heart rate variability represents a mathematical correlate for adapting the heart rate to changing requirements in the human organism and is regarded as an expression of neuro-vegetative regulatory capacity. The RR intervals are e.g. measured in milliseconds [ms] and fluctuate i.d.R. between about 700 and 1200 ms.
The determination of heart rate variability (HRV) and associated physiological characteristics, i.a. Heart rate, QPA, VQ, SNNNRR, is well known in the art and makes it possible to detect subtle changes in the regulatory system of the human body. The autonomic nervous system, also called vegetativum, regulates, among other things: heart activity, blood pressure, distribution of bloodstreams, respiratory depth, respiratory rate, thermoregulation, glandular secretion and gastric and intestinal motility. It is divided into two subsystems, the sympathetic and the parasympathetic nervous system.
The heart rate or heart rate is determined on the basis of time indices of the detected R-R intervals. To describe the tone of individual regions of the autonomic nervous system, a spectral analysis is performed (FIGS. 3-6). The measurement data of the heart rate are in this case transmitted by means of mathematical methods from the time domain to the frequency domain and represented as a power spectrum. Thus, various other frequencies present in the heart rate can be made visible, which are empirically assigned to specific body rhythms such as respiration or blood pressure. In principle, all HRV sensor and software systems known in the art or to be developed in the future may be used for the analysis method according to the invention which is described below. Statistical processing of received data records takes place e.g. with proprietary evaluation routines such as MatLab® and SPSS®. For the calculation of a frequency analysis, the heartbeat sequence is divided into equidistant sections and brought by a Fourier transformation from the time domain to the frequency domain. This transformation decomposes the total signal into individual sinusoids and reproduces the magnitude of the individual frequency components. The frequency bandwidth of substantially 0 to 0.5 Hz is calculated according to the order of magnitude of the frequency portions for specific time periods. In particular, the following frequency ranges are distinguished: 0.04 to 0.15 Hz: Low Frequency (LF) with correspondence: predominantly sympathetic activity, to a lesser extent also vagus activity, in particular allocation of the blood pressure and circulation rhythm; 0.15 to 0.40 Hz: High Frequency (HF) with correspondence: parasympatheticusA / agusactivity; Allocation, in particular of respiratory functions, corresponds in particular to the parasympathically determined oscillation portion of the respiratory sinus arrhythmia (RSA) and thus to the respiratory-synchronous cardiac frequency fluctuation;
The power over the entire frequency range of 0.0033 - 0.5 Hz is called TOT (total frequency). The TOT indicates the size of the total area within all frequency ranges and is considered the measure of the influence of the vegetative on the organism. Within the respective frequency ranges, the power is determined and converted using the natural logarithm. Further, as special frequency ranges, ULF (ultralow frequency; < 0.003 Hz) and VLF (very low frequency, substantially 0.003-0.04 Hz) can be cited.
The activities of respective frequency regions or vegetative organ functions over a given period of time are stored in an HRV frequency information record, which in practice can be represented graphically as a spectrogram (see Fig. 10). The spectrogram is a clear representation of the complex rhythm information contained in the heart rate or heart rate variability. In chronobiology, it is considered an image of human proper time and is also called an "AutoChrone image" ("autos" = self, own; "chronos" = time), where the information is in three dimensions (here: abscissa = time, ordinate = Frequency, color = amplitude). The amplitude indicates the energy density of respective frequency ranges determined from the HRV frequency analysis and is preferably stored as a dimensionless numerical value. The scale amplitude can be chosen arbitrarily and e.g. from 0 to 1 or from 0 to 50, where direction zero-tending values denote a very low amplitude strength and directions 1 and 50, respectively, a very high amplitude strength.
Each (in Fig. 10: vertical) line of the spectrogram is the result of the frequency analysis of a short section of a time series, e.g. a heartbeat episode. The amplitude of the respective rhythms is coded in color. A small amplitude is e.g. shown in blue, a higher in white and yellow, a very high in red (see also Fig.6). It should be noted that the original color graphics according to the figures were each resolved into pure black and white puncturing and therefore the color coding of the amplitude intensities is only poorly recognizable.
In particular, the spectrogram allows a statement about the sleep architecture, whereby a good sleep is cyclical and quiet sleep phases clearly differ from REM sleep phases (dream sleep). A bad sleep appears fragmented and vegetatively restless. The vegetative balance (tension recovery, see also the following defined vegetative function parameter VQ) in good sleep is vagotter than in bad sleep. Using a frequency analysis, the heartbeat and the spectrogram can be used to calculate autonomic rhythms and sleep stages. These meta-rhythms develop from the circadian rhythm, which is characterized by the change of day and night. The day and night events then follow the BRAG, a cycle of rest and activity cycles in a 2-h rhythm associated with the circadian rhythm. Activity phases (90 - 120 min.) And regeneration phases (5 - 20 min.) Alternate during the day, which continue in the nocturnal sleep architecture as longer deep sleep and shorter REM phases. The perfusion of peripheral tissues reveals a minute rhythm (= 0.017 Hz), the rhythm of the blood pressure follows a 10 second oscillation (= 0.1 Hz), people who are tensing or sleeping show the rhythm of breathing (= 0.2-0, 3Hz) in the heartbeat.
Other important parameters of the autonomic nervous system or vegetative function parameters are listed below:
Heart rate: This characteristic already mentioned above denotes the absolute interval duration between two R-waves or their differences. This results in the representation over the time axis to the heart rate gradient. The heart rate is thus the number of heart beats - represented by the R-waves - per unit of time, usually per minute. In the consideration of the HRV, however, it should be taken into account that at a heart rate of e.g. 60 beats per minute can very well change the heart rate within that minute. Thus, the heart rate should be considered as an average value within one minute. However, if one does not consider the absolute number of R-waves in the time domain but the distances from R-wave to R-wave (see FIG. 1), one arrives at a time measure of the time differences of the RR-distances in milliseconds (ms). , Analysis in the time domain of absolute RR intervals results in heart rate variability (HRV) - that is, the variability of heartbeats. LF / HF (vegetative quotient, VQ): The quotient of the two frequency ranges LF and HF contained in the aforementioned HRV frequency information data set reflects the current vegetative activation level of the organism. Higher VQ values show an active, performance-oriented body setting, low VQ values a recovery-oriented one. The VQ is represented over time and represented as the ratio of activation of sympathetic (preferably red) and parasympathetic (preferably blue). A VQ of 1: 1 thus means an equal occurrence of sympathetic and vagal rhythms. Falling below this, with a line-marked threshold, can be equated with the onset of a good recovery; If the test subject also remains above the ratio VQ = 1: 1 in the sleep phase, the recovery is disturbed. SDNNRR (standard deviation of normal-to-normal intervals): The standard deviation above i.d.R. 5-minute artifact-adjusted RR interval series is a measure of overall variability over all frequency ranges. Long-term studies have shown that overall variability is a measure of vitality and that people with lower SDNNRR live much shorter. Statistically, the SDNNRR is the measure of spread around the mean of the interval duration (or its difference) within a predefined period of an HRV analysis. This will show the amount of variability of all RR distances - and hence heart rate variability - over a specified period of time.
The standard deviation is generally calculated:
Respiratory Sinus Arrhythmia (RSA): the respiratory change in heart rate. The RSA calculates to RSA = median (| HRi-HRt-1 |), where the median represents the central value of the sorted number sequence within a time period. The RSA represents the dependence between pulse and
Breath and shows in the spectrogram as an oscillation between 0.2 and 0.3 Hz, which occurs during sleep and relaxation phases. This phenomenon is caused by the respiratory-synchronous fluctuation of the heart rate, which increases with inhalation (sympathetic) and decreases with exhalation (vagus). The function of RSA is to maximize gas exchange through an optimal interaction of perfusion and ventilation with each breath. The RSA is converted to logRSA using decadic logarithm. logRSA: The median of the absolute differences of consecutive heart rate values, like the HF, primarily refers to the rapid, breath-induced changes, but without drawing a strict limit at a particular frequency. Respiratory sinus arrhythmia (RSA) is thus the high-frequency variability of heart rate, which reflects the rate of cardiac rhythm modulation by respiration. It is also a measure of the tone of the vagus activity. QPA: The pulse-to-breath ratio indicates how often the heart beats during a breath (ratio of heartbeats to a respiratory cycle). During the night and at rest it could be observed that in healthy persons a QPA setting comes to a ratio of about 4: 1, regardless of the quotient under ergotrophic conditions during the day, which may be between 2: 1 and 22: 1.
In the course of a method according to the invention for analyzing the interaction of electromagnetic radio-frequency immissions (HFI) with vegetative regulatory mechanisms of a test subject, the following method steps are fundamentally provided: Determining the course of an intensity of radio-frequency immissions (HFI) measured as power flux density with respect to the test subject over a defined one Period, - determining the variance in the time interval of consecutive heartbeats detected by electrocardiogram (ECG) measurement and the progression of the hereditary variability (HRV) associated vegetative function parameters (HRV-FP) of the test subject over a defined period of time determined therefrom, - wherein the HFI and the HRV-FP profiles are at least partially synchronized and calibrated to at least one common reference time interval T. In practice, the measured values acquired by the ECG measuring device or HRV-FP evaluation device and the HFI measuring device are measured. Data device provided in different recording frequency. In order to allow a computational comparability of the ECG or HRV-FP data and the HFI data, interpolation and / or extrapolation and / or selection for both (HRV-FP and HFI) data channels to be compared in each case an equal number of measurement or Thus, in the evaluation practice, temporal synchronization means that measured values or data of both (HFI and HRV-FP) signal channels are chronologically indexed and quantitatively coordinated with one another in the considered reference time interval. Performing an analysis of the HFI history in which primary significances 1 -n are detected, which are temporary, i. within the reference time interval T ^ n, time-limited deviations of the HFI curve from defined static or dynamic mean values or HFI reference data, in particular significant increases and / or decreases in the HFI curve compared to immediately preceding and / or subsequent sections of the HFI curve; In any case, the primary significances, as well as the secondary significances described below, in any case last shorter than the entire reference time interval Τ ·, _η. Performing an analysis of the HRV-FP history, Tn which secondary significances T-n ', Le.temporary, i. within the reference time interval Ti.n time-limited deviations of the HRV-FP curve of HRV-FP reference data or rule progressions, which is derived from empirical data and / or dynamically generated, i. are derived from current measurement / calculation data or HRV-FP progressions, in particular significant increases and / or decreases in HRV-FP progression with respect to temporally immediately preceding and / or succeeding sections of the HRV-FP run; upon detection of a defined number of coincident, i. in the HFI course, and in the HRV-FP course, a positive evaluation of the presence of a bumanbiologically relevant influencing of the vegetative regulatory mechanisms of the test subject takes place substantially simultaneously occurring primary and secondary significances 1-n, 1'-n '.
Changes or correlations in the HRV-FP course and / or in the HFI course are immediately and immediately, i. recorded in real time. The measurement of the HFI intensity or power density is carried out, for example, with a Maschek ESM 140 dosimeter, but can also be carried out with other suitable HFI measuring devices and spectrum analyzers. High-frequency emission measurements can in principle be made with broadband meters or with frequency-selective equipment. Broadband meters are easy to use and allow quick overview measurements. However, because of their low dynamic range and, above all, the lack of frequency selectivity, they are not suitable for precisely determining from the total emission spectrum that part of specific high frequency emitters such as e.g. Mobile base stations goes out. In addition, no correct extrapolation of the measured power flux density instantaneous value to the immission at maximum system utilization is possible.
Within the evaluation, the originally measured power or voltage values are converted into the power flux density or field strength and extrapolated to the HFI at maximum system load. This is possible with GSM base stations by identifying the immission of the control channels BCCH from the measured spectrum and multiplying by a factor which is determined by the number of maximum available frequency channels. It is customary to specify the power flux density in units W / m2 or W / m2 or in mW / m2 or pW / m2 (1 watt = 1 W, 1 milliwatt = 1 mW = 0.001 W, 1 microwatt = 1pW = 0.000 001 W).
Whenever we talk about power flux densities (LFDs), we refer to HFI intensities at specific times. The terms HFI and LFD are thus to be understood synonymously. Of course, the HFi waveform diagrams scaled in the following figures 7, 9, 12 in unit mW / m 2 could also be presented in alternative units for measuring the power flux densities caused by external high frequency sources. The measurement of the power density is continuous, individual HFI Measured values are provided with time indices and logged chronologically. The recording intervals of the HFI meter are e.g. 0.5 seconds.
The HFI meter is placed directly on the body, e.g. attached to the arm or in the immediate vicinity of the body of the test subject. It preferably has a frequency-selective measurement of mobile radio bands corresponding to the state of the art or of other frequency ranges used by HF wireless technologies, in particular GSM-900, GSM-1800, DECT, UMTS, LTE, WLAN and the like. The measurement data acquired by the HFI meter can be output as a graphical or tabular representation. The recording of the HFI meter is initiated by the test subject by operating a switch or other input method. Calibration and synchronization with the recording data sets of an HRV-FP measuring device also provided for the observation of the same reference time interval, but time-delayed activated and deactivated, is carried out by a method-conforming evaluation algorithm based on respectively assigned time indices of the HFI and HRV-FP data sets or event histories.
Detected HFI frequency ranges or mobile radio bands can be selectively included or hidden in a subsequent evaluation. First of all, the use of those HFI courses or mobile bands in which the highest power flux densities or the most primary significances were determined during the investigation period is recommended. In addition to the standard power density average RMS (6 minute ICNIRP averaging) average, peak function is also possible to represent the actual peak HFI loads over time. For each measurement sequence, a unique identity code is defined which later enables the identification and timing of the performed HFI measurements. For graphical evaluation, selected frequency bands or HFI waveforms can be exported as an image file (e.g., using the Visual ESM-140 visualization software).
In particular, the following characteristics are detected as primary significances 1-n in the HFI curve: maximum values or LFD peaks, which in the case of a graphical representation can be seen as approximately needle-shaped or jagged amplitude excursions, approximately step-shaped course segments, in this case preferably the initial and / or end sections of a respective approximately horizontal or a constant LFD indexing run section of the staircase form and / or the same approximately horizontal course section itself, - several, preferably more than three, preferably within a time interval of maximally 10, 30 or 60 minutes successive increases and decreases of the course, Exceeding an LFD of more than 0.1 mW / m2, preferably of more than 0.05 mW / m2, more preferably of more than 0.01 mW / m2, - approximately constant LFD or in the case of a graph substantially parallel Timing sections d it HFI course, - subsidence and / or increases in the HFI course or the LFD, which preferably have a difference of more than 10%, particularly preferably of more than 30% compared to immediately preceding and / or following HFI-course sections; as immediately preceding 1 may be a period of several, e.g. 1-30 minutes, preferably 0-60 seconds, more preferably 0-10 seconds.
As primary significances, not necessarily only HFI / LFD peaks or the occurrence of the aforementioned characteristics need to be defined. On the other hand, a phase of the HFI process in which the LFD is particularly low and thus less burdensome can also be a factor. The detection criteria for the detection of primary as well as secondary significances can thus also be defined negatively, so that a significant improvement of vegetative functional parameters (preferably in the case of simultaneously detected change or immission reduction in the HFI curve) is also detected as a secondary significance, which indicates a human biologically relevant influence of the test subject ,
The mathematical calculation is based on an HRV diagram sequence illustrated in FIGS. 16-19 explained below, from which exemplary comparison points were selected over all time intervals, which each coincided with a primary significance in the HFI curve and subsequently compared with HRV-FP progressions determined at the same time were.
For analysis evaluation, detected primary significances 1-n may be differently weighted according to their specific characteristics and / or their occurrence in respective HFi frequency ranges. For example, a relatively small, spontaneous increase in HFI progression after a phase of relatively low or equal HFI stress may result in a more significant effect on concurrent HRV-FP progression than a sharp increase in HFI progression in already mediated LFD Peak times T1H1 and HFI loads, respectively. Thus, as well as the height of the LFD peaks, the size of the current distances between significant changes in the LFD or between primary significances 1-n can also be used as an assessment criterion, with events or significances with longer intervals assigned a greater burden prevalence than those with smaller time intervals.
Since mobile wireless transmission standards each have specific clock rates or pulse frequencies, it may also be the case that the human organism is limited to individual HFI
Frequency bands or mobile bands react more sensitively than others. In analysis practice, it has been shown that HF1-conditioned physiological reactions of the test subject in the form of inventively selected secondary significances 1'-n 'individually precipitate, and therefore only a limited significance is given to a measurement of mere HFI power flux densities customary in the interior space sector. For example, it may be the case that the test subject responds to an LFD of 0.02 mW / m2 in the WLAN frequency band with equally high physiological load as to an LFD of 0.5 mW / m2 in the GSM frequency band. An analysis of the HRV-FP curves proposed in accordance with the invention in conjunction with HRV-FP reference data and HFI courses avoids this evaluation problem and makes the respectively given, HF1-related exposure of the vegetative organism of the test subject directly visible. Exemplary assessment factors or prevalence factors for evaluating primary and secondary significances detected in HFI traces and / or in HRV-FP traces may be stored in algorithmic form on a memory device and used as automated processor routines.
To determine the ECG data or the HRV-FP courses, a suitable ECG recorder with associated HRV evaluation unit is used. Due to its compactness and portability, the HRV measuring device "HeartMan" from the company HeartBalance ® has proved its worth in the present field of application due to simultaneous measuring precision. The HRV meter comprises three self-adhesive electrodes E1, E2, E3 and, due to its small size, can also be attached to the body of the test subject during sleep, e.g. be glued by plaster strip.
As HRV-FP, all can be termed "time domain measures". or "Frequency Domain Measures" known or modified parameters suitable for constituting a time course. It is understood that it is also possible to use the parameters or HRV-FP standardized in the course of HRV measuring methods, e.g. by statistical operations or by a variation of the number or recording frequency of each considered measurement or calculation parameters to be adapted. In particular, if in the present context of SDNN is mentioned, this statistical characteristic could also by equivalent or adapted Statististische Kenngrößenbzw. "Time domain measures" are replaced as e.g. SDNN (SDNNIDX / ASDNN), AVNN (Average of all NN.), SDNNIDX / ASDNN (Mean of the standard deviation in all 5-minute segments of a 24-h recording) intervals), rMSSD (Square root of thematic of the squares of the differences between adjacent NN intervals), pNNx such as pNN50 (Percentage of difference between adjacent NN intervals> 50 ms)
The data recorded by the HFI meter and the HRV meter (electromagnetic power flux densities and RR intervals) are recorded with different recording frequencies depending on the device. The HFI meter or dosimeter measures and stores the value of the current power flux densities approx. Every 6-7 Seconds, while the heartbeat or the resulting RR intervals are recorded in a much higher recording frequency. In order to make the data collected by both measuring devices (if applicable, more than two HFI and / or HRV measuring devices can be compared with each other), these must be converted to the same time units so that the same number of measured values is found in each of the two signal channels. For this purpose, the reference time interval for both signal channels or for both (HFI and HRV-FP) courses is subdivided into a plurality of time blocks or sub-time intervals. The duration of a sub-time interval can be configured variably as input parameter depending on the requested accuracy and characteristic of the analysis. It is recommended to create sub-time intervals of 1-10 seconds, preferably 3-8 seconds, to preserve as many RR interval readings as possible provided by the ECG meter. To achieve additional accuracy, overlapping sub-time intervals may be calculated become. Subsequent to the sub-time intervals, various operations (e.g., selection / detection of maximum, minimum, median, average or delta values) may be applied to detect the most interesting information according to the significance of the signal to be considered. In the case of the HFI curve, for example, the maximum value of the sub-time intervals is selected in each case. In order to ensure a high resolution and detection accuracy, the sub-time intervals in the present embodiment are 7 seconds each.
As secondary significances 1'-n 'in the HRV-FP curve, it is possible to detect all events which, when plotted, have a slope k> 1 or a pitch angle α of more than 45 °, preferably of more than 70 °, particularly preferably of have approximately 90 ° with respect to the time axis (slope k- * °°) or one or more approximately needle-shaped elevations of the HRV-FP curve. Analogous to the assessment of the HFI curve, suitable detection criteria can also be defined and weighted differently with regard to the HRV-FP curve. As a detection criterion for secondary significances, e.g. an exceeding of the average HRV-FP curve or of the HRV-FP curve section, which is analyzed in each sub-time interval, directly and / or subsequently HRV-FP curve section should be set by at least 20%, 30%, 50% or 100%. Since the autonomic functions of humans are governed by a variety of autonomic rhythms, any deviating arrhythmia or distortion of a HRV-FP history can be defined as a detection criterion for a secondary significance, and in particular all The events illustrated in the following analysis examples.
Numerous modifications and specifications are possible in order to optimize the validity of the analysis method according to the invention or to minimize the required measurement duration. In a specific method variant, for example, a targeted initiation of a radio-frequency immission acting on the test subject can take place by means of a specially provided test radio-frequency source and thus one or more temporally determined primary significances 1-n can be generated. A time-averaged HRV-FP curve is checked for coincident secondary significances 1 '-n'. The radio frequency emissions generated by the test radio frequency source may preferably have realistic power flow densities typically provided by transmitters according to conventional radio frequency standards such as GSM, UMTS, LTE, WLAN, Bluetoooth / IEEE 802, DECT and the like. in general, especially in the domestic area, correspond or simulate those Hochfrequenzquetten.
Also, the test subject may during one or more time intervals by a preferably surrounding on all sides of the test subject shielding device such. be shielded by electrically conductive metal fabric of external high-frequency emissions of the environment. The shielding device is either opened for a defined period of time or, with the shielding device closed, at the time interval of the shielding, a test radio frequency source located within the shielding device, e.g. a mobile phone activated one or more times for a defined period of time, wherein currently determined HRV-FP history of the test subject coincides with the secondary significance (s) generated by the test radio frequency source, or during the time period of the test subject Opening the shielding occurring secondary significances is checked.
It has e.g. proven, the test subject 5-10 min. in a first phase, with the shielding device still open, to be exposed to HFI, in an immediately following second phase, the shielding device for 5-10 min. close and in a subsequent third phase the Abschirmvorrichtung for 5-10 min. open again (the specified duration may also vary). Here, in one or more of the stated phases, activation of the test radio frequency source, e.g. made of a mobile phone. Such activation is preferably done without the knowledge of the test subject, either automated or initiated by a third person. The test subject wears a blindfold during testing to prevent falsification of the test result by visual stimuli.
In a further development of the analysis method according to the invention, provision may be made for a vibro-acoustic relaxation program to be carried out before the test object is tested, preferably during the first phase described above, in which the test subject is exposed simultaneously to two different but correlated actions: on the one hand, a vibration sequence generated by a lying with mechanical-electrical actuating means vibratable lying device, on the other hand, an acoustic sequence. The vibration and / or the acoustic sequence is in this case modulated on the frequency corresponding to the BRAC (basic rest-activity cycle). Previous experiments have shown that using a particular BRAC rhythm-based vibro-acoustic relaxation program provides a particularly high level of detection accuracy of HF1-induced secondary significances. As demonstrated in a pertinent publication, emotional-functional enhances after performing the vibro-acoustic sequence Stress parameters / EFSP significant (Matissek, 2011: "Follow - up of the mood change of 20 subjects using acoustic rather than light - induced stimulation of the base - resting activity - cycle / BRAC"). The use of such relaxation program, which is already known from the prior art, has therefore proven to be particularly advantageous in combination with the analysis methods proposed according to the invention and focused on HFI loadings. By exposing the test subject to two different but correlating sensory stimuli, any additional external stimuli which may affect the detection of the HRV-associated autonomic function parameters may be largely eliminated.
According to a further variant of the method, depending on the determined frequency and / or standard deviation (size) of the HFI-induced secondary significances, e.g. in percent, gradual assessment of the likelihood of intolerance to radio frequency immissions. Also, a gradual evaluation of the vegetative regulatory capability of the test subject against HFI-induced secondary significances can be made by comparison with HRV-FP comparative data from external test series. The analysis according to the invention or the comparison of the event histories of the HFI and / or HRV-FP profiles over a selected reference time interval can each be automated, semiautomatized or also manual. An illustrative comparison results, in particular, in the case of a graphical overlay of HFI and HRV-FP progressions, which were at least partially synchronized and calibrated to at least one common reference time interval T 1.
The subsequent diagrams exported as digital graphics of the comparative evaluation always start on the left side with the same start time of the comparison period. Subsequently, all the other graphics, in particular HRV-FP progress diagrams and activity protocols, can be positioned exactly one below the other. In order to allow the client or the test subject a fast and plausible overview of the correlations detected in the evaluation or the correlation between primary and secondary significances 1-n, 1'-n ', the progress diagrams can each have a preferably orthogonal to the time axis Thus, in the present embodiments, vertical grid lines may be provided. Both equidistant raster lines and specific raster lines can be provided, which can be used at times relevant events or when primary and secondary significances 1-n, 1'-n 'occur in the HFI and / or. HRV-FP curves were set and can also extend over several diagrams according to Figure 7. In this way, the immediate correlation of a trigger in the HFI curve with corresponding effects in the HRV-FP curve is illustrated. The HRV meter and / or the HFI meter may also be provided with a GPS locator function to provide a verified record of the local relationship of the analysis result.
Alternatively, to establish an immediate correlation of a trigger in the HFI course with reactions in the HRV-FP, the determination of such a correlation may also be indirect:
In a special process variant, an automated comparison of the currently ascertained HRV-FP curve with characteristic, HF1-induced secondary significances or with HRV-FP reference data or control curves is carried out in each case from the current experimental setup or from at least one of these preceding test arrangements with the current test subject or third test subjects have been determined and are now held in the form of digital or analog data, preferably in the form of graphical histories, alphanumeric or algorithmic data, functions or function graphs on a memory device. Upon detection of a defined, at least partially present agreement of the currently determined HRV-FP curve with characteristic, HFl-induced secondary significances or deviation of the currently ascertained HRV-FP curve from the HRV-FP reference data, a positive assessment of the presence of a human biologically relevant curve is again obtained Influencing the vegetative Regulati¬onsmechanismen made the test subject. All arrhythmias or secondary significances Tn 'documented on the basis of subsequent embodiments could therefore also be detected on the basis of aforementioned recognition algorithms or in comparison with characteristic, HF1-induced secondary significances or HRV-FP reference data obtained at an earlier time in the course of HFI correlation analysis ,
The diagram according to FIG. 8 shows a heart rate profile of the test subject over a period of eight hours. In the history diagram (FIG. 7) arranged above it, the simultaneously measured HFI profile in the mobile radio network GSM900 is considered. The graphs were calibrated to the sleep phase of the test subject, starting with mark A and ending with mark D (sleep end).
The diagram according to FIG. 8 as well as all HRV-FP profiles described below can each also be provided with activity protocols, so that spontaneous events such as e.g. a possible awakening and interruption of the sleep phase may be associated with corresponding abnormalities in the HRV-FP courses, or in this case detected secondary significances are eliminated as non-HFI-induced from the evaluation. According to FIG. 8, for example, sleep abruptly stopped at mark B at about 4:00, but resumed shortly thereafter. In the morning, the test subject wakes up at mark C (about 7 o'clock), but sets the sleep again to mark Dfort, but then wakes up again and gets up in the episode (7.30 o'clock).
Excursus on heart rate variability (HRV): The dark main line of the heart rate shown in the heart rate diagram according to FIG. 8 corresponds to the average of the heartbeats, represented by the R waves, in the time interval of one minute. However, there is also a metrological determination of Intervalldau¬er the individual heartbeats relative to each other, which in diagram form as fine, in transverse to Zeitachsebzw. vertical direction above and below the dark average curve. Lines running across the average curve each indicate a shortening of the heartbeat interval between two heart beats, and lines extending below the average curve each extend the heartbeat interval. A visualization of the heartbeat intervals shown in chronological order, as shown in Fig. 8, appears in an approximately tubular meandering passage. The average and interval-duration curves originally available in color or in light-dark color were resolved into pure black-and-white puncturing for the present documentation. The resultant pixelated representation thus does not correspond to the original progression graphic (this note also refers to all the figures described below which are originally color-coded in the original).
In principle, one can describe three characteristic courses in the heart rate diagram: "Machine cycle": The heartbeats follow an exact regularity and thus show a limited adaptability of the organism. There are only very short vertical thin lines above and below the average curve (lower HRV). - Normal Variability: The variation of the interval duration from heartbeat to heartbeat shows theregulatory capacity of the organism. The farther the representational lines extend above and below the mean average curve, the faster and more flexible the heart adapts to internal and external influences (high HRV). Significant arrhythmias / extrasystoles (double strikes): The interval between bursts is extremely shortened, which is shown in the graphical evaluation in a vertical long line above the average curve and thus in the form of a significant deviation of the adjacent variability. The longer the objective line, the more it deviates from the average interval duration.
In the course of the method according to the invention, it is provided that the time course of the heart rate is analyzed and in this case the occurrence of short-term arrhythmias, i. Irregularities are detected as secondary significance 1'-n 'to each of the temporally preceding and / or subsequent heart rate progressions. The detection accuracy can be determined parametrically or according to a gradual determinable sensitivity factor, so that arrhythmias of greater severity can be recognized as secondary significance, such as weak, masked, distorted or partially compensated arrhythmias. For this purpose, the heart rate course as well as other HRV-FP courses can be viewed under any macro or micro scale or scaling. The same applies to the detection of irregularities in other HRV-FP described in the present application. As a particularly significant arrhythmia in heart rate progression, extrasystole is preferably detected.
To Eventislorlorie In the courses according to Figure 7 and Figure 8: The ZeltintervaHTI shows a very low. negligible LFD in the HFI diagram according to FIG. 7 and likewise no extrasystoles in the heart rate diagram according to FIG. In the time interval T2, a sudden increase in LFD occurs, which is characterized by a chaotic sequence of short and substantially constant LFD (staircase shape) and small HFI peaks. It has been shown that such a scenario exerts the strongest stimulus on the organism, which reacts immediately with a massive increase in the heart rate. The intervals T3 to T5 show the typical course of consistent LDF (staircase form). This always begins and ends with an LFD peak (primary significance 2-5), which also regularly shows the simultaneous occurrence of an extrasystole and thus a secondary significance 2'-5 'in the heart rate diagram. The organism of the test subject is less activated by these LFD tips and the subsequent concomitant LFD in that it compensates for extrasystoles immediately in normal variability (each time nearly equal lengthening of the heartbeat interval) and the lower level of variability of LFD in overall slightly increased variability in HRV .This also produces corresponding stair patterns in the course of the heart rate (Fig.8).
Time interval T6 shows substantially similar events as the time intervals T3 to T5 with respect to HFI and LFD characteristics, respectively. However, the time span between the constant LFD and the LFD tips is about half as long as before, with the LFD peaks slightly larger. This also leads to the increase in extrasystoles in exact time sequence and an increase in the heart rate (Figure 8). Time interval T7 is again characterized by a relatively low LFD but with 2-3 LFD peaks. Overall, the heart rate lowers again and offsets the externally triggered extrasystoles (LFD peaks or primary signatures 14, 15) by the increased normal variability. In time interval T8, there is an extremely short change between LFD plateaus and LFD tips, which have also increased in height. The heart rate in this case corresponds to the LFD, increases sharply as a result of the short LFD change, with a corresponding increase in extrasystoles, and immediately falls again as the LFD decreases, in order to rise again at the subsequent LFD peaks (compare primary significance 14-22 in FIG of these induced secondary significances 14'-22 'in Fig. 8). Eventually, this LFD trajectory results in a sharp, sudden increase in heart rate associated with the same-time awakening (mark B) that will repeat in similar scenarios, as illustrated below by mark C and D.
Time interval T9 again shows a very low and constant LFD (FIG. 7), the heart rate also decreases in the sequence (FIG. 8). The long vertical lines below the dark heart rate average curve indicate the active regulatory activity of the organism, which again leads to a continuous lowering of the autonomic stress level (corresponding courses are further explained by further parameters such as QPA and SDNNRR / logRSA). Time interval T10, like time interval T2, shows a very chaotic sequence of very short, consistent LFD and high LFD peaks (see primary significances 25-36 in Fig. 7). The heart rate variability is very high in the follow-up revascularization with correspondingly frequent occurrence of extrasystoles (see secondary significances 25'-36 'in Figure 8).
The HFI curve or the LFD in the time interval T11 initially exhibits a similar characteristic as in the time interval T10, but with increasing height of the LFD peaks. At marker C, a simultaneous occurrence of a very high LFD peak or primary significance 39 and a temporally exact extra systole 39 'is detected. At this time, the test subject wakes up briefly, falling asleep again, shortly thereafter awakening at the same time again with very high LFD peaks or primary significances 40, 41 and, respectively, with temporally exactly matching extrasystoles or secondary significances 40 ', 4T. It is followed by another fall asleep and finally final waking, which in turn coincides with a simultaneously detected LFD peak 42 and associated extra-systole 42 ', respectively. The coincident sequence of primary significances 1-42 in Figure 7 and secondary significances 1'-42 'in Figure 8 shows that a randomness in the correlation between HFI loads and heart rate progression or occurrence of extrasystoles is excluded. Thus, in the previous comparison in a reference time interval of 8 hours, there is a correspondence between primary and secondary significances of approximately 100%. The externally induced activation of the heart rate leads to an increased burden on the heart, which is also expressed in an increased number of strokes. This can be calculated exactly from the heart rate, whose level in the time interval Ti (almost no HFI / LFD) was about 53 beats per minute. Since the heart rate increases only slightly and in the short term during sleep in the REM phases (see also, inter alia, explanations regarding sleep architecture and RSA), there is an additional burden on the heart of approximately 6.878 strokes in the reference time interval.
Under the given conditions, the organism can no longer follow its own vegetative-controlled rhythm but is forced to compensate for the events of external, HF1-related influences permanently, which leads to a disturbance of the body's own regulatory processes. Such a disorder of the regulatory processes may lead to the EMF symptoms mentioned above, such as diminished performance, sleep disorders, poorer regenerative ability, mood disorders, anxiety and depression. While an assignment of such nonspecific disorders to respective EMF exposure had previously been tested by environmental medicine outpatient clinics in laborious differential diagnostic procedures, blood findings and laboratory analyzes, it is now possible to make such an assignment in a time and cost-saving manner.
Additionally or alternatively, an influencing of the test subject described above by HFI can also be detected on the basis of further vegetative functional parameters (HRV-FP), e.g. by an analysis of the already initially defined pulse-to-breath ratio (QPA) according to FIG. 12, a plurality of primary significances 1-29 were again detected, which correlate with substantially contemporaneous secondary significances T-29 'in the form of QFA maxima in the QPA performance diagram of FIG. The history of the event in detail:
Marks A, B, C, D again correspond to the events: going to bed at 23:57, interim waking up and falling asleep at 4:00, waking up at 7:00 and finally getting up at 7:30. The time interval T1 shows almost no LFD in the HFI curve according to FIG. 12, the QPA according to FIG. 13 being approximately 4: 1, which is considered to be the optimal ratio and the indication variable for a normally functioning regulation. In the QPA diagram of Fig. 13, therefore, the level 4: 1 is marked with a horizontal broken line.
In the time interval T2, a sudden increase in LFD begins, which is characterized by a chaotic sequence of short and consistent LFD and a small LFD peak (detected as primary significances 1-4). The organism responds to this scenario with a significant increase in QPA (see especially secondary significances 2 ', 3', 4 'in Figure 13). The QPA curve thus jumps from the optimum ratio of 4: 1 to 7: 1 within a few minutes and then drops to 6: 1 (secondary significance 4 '). The intervals T3-T5 show the typical course of consistent LDF (staircase shape).
This always begins and ends with an LFD peak, which also regularly shows up in the simultaneous occurrence of a temporary peak in the QPA course. If the LFD falls, then the QPAab begins to sink as well. In the time interval T6, T8, T10 and T11, there are further significant correlations. Each LFD peak in Fig. 12 causes a QPA peak and secondary significance in Fig. 13, respectively, with the QPA decreasing immediately after reduction of the LFD.
The magnitude of the QPA peaks, or secondary significances, is variably related to the strength of the LFD peaks and is particularly enhanced by extremely short changes in LFD in the second or second cycle. For example, In the time interval T8, there is an extremely short change in the HFI curve or in several short-time successive LFD tips, so that the QPA curve increases 14 'higher in the case of the secondary significance (8: 1) than in the case of LFD tip-like strength (eg secondary significances 6 ', 8', which "only" lead to an increase of the QPA to the ratio 5: 1.) In the time interval T11, this relationship is again clearly visible at the highest LFD peaks (primary significances 27-29) 9: 1 secondary significances (secondary significances 27'-29 '), again showing a significant match between LFD peaks in the HFI curve (Figure 12) and QPA peaks in the HRV-FP curve (Fig In this comparison, in a reference time interval of 8 hours, there is a match between primary and secondary significances of approximately 100%.
The diagrams according to FIGS. 14 and 15 again illustrate, analogous to the preceding FIGS. 12 and 13, an analysis of the QPA profile, but now with regard to lowering of the QPA (the observed time period T1-T11 or the HFI and QPA profiles are with those as shown in Figs. 12 and 13, however, were provided with new, topic-specific markings).
This is followed by a temporal comparison of all time intervals with low LFD with the QPA course. As mentioned above, in a healthy organism, a QPA of 4: 1 sets in during sleep. If external events such as LFD peaks trigger an activation, the organism immediately afterwards strives for the optimal balance ratio (synchronization). The time interval T1 indicates the starting situation with a QPA of 4: 1 without HFl-related activation and thus still normal functioning regulation. Already a consistent LFD does not produce any or only a small stimulus and the QPA starts to decrease in the sequence. In time intervals T3 and T4, a period of 15-20 minutes with LFD remaining constant results in the return of the QPA profile to the ratio of approximately 4: 1. This is repeated in the time intervals T7, T9 and T11. In the time interval T8 and T10, a chaotic sequence of short-steady LFD levels and LFD peaks results in a 5: 1 permanent increase in QPA, which is no longer below.
As the above example shows, a very unstable equilibrium ratio between organ groups such as the heart and the breathing synchronized with it in the sleep phase (but also other organ groups or vegetative control circuits are affected) is impaired by very low LFD, which is essential for the regeneration in the deep sleep phases Meaning is. Numerous studies on HRV have shown that even with very different daily outcomes of QPA, there is always a normalizing effect on the 4: 1 ratio at night, with the vegetative mechanisms during sleep as a system of highest order versus chaotic
Apply to daily routines. As it turns out, this order, visible in the integer coupling, is completely lost by an external, HF1-related activation and a subsequent, chaotic course of the OPA. The optimal integer ratio of the OPA of 4: 1 is shown only at more time intervals with almost no LFD or HFI load (see, e.g., T1, T4, T9 with significances 1 / T, 3/3 ', 7/7'). Thus, in the present consideration, not the LFD peaks and temporally corresponding QPA peaks were evaluated as primary and secondary significances, respectively, but approximately plateau-shaped waveforms, which temporally correlate in the HFI course and in the QPA course. Such evaluation may be complementary or alternative to the previously described detection of correlating peaks in the HFI and QPA.
Incidentally, the 4: 1 QPA regulatory structure occurring during sleep can be determined not only between heartbeat and respiration but also with regard to blood pressure and peripheral tissue perfusion rhythm. An analysis according to the invention with regard to correlating significances could thus also be carried out analogously by consideration of functional parameters or HRV frequency ranges based on blood pressure or peripheral tissue perfusion rhythm.
As a further HRV-associated vegetative function parameter, the time course of the SDNNRR (standard deviation of normal-to-normal intervals) of the test subject can be used. An SDNNRR waveform shown in FIG. 21 is a representation of a time-spread statistical spread about the mean of the heartbeat interval duration or its differences. This function parameter is suitable for representing the change in substantially all frequency ranges of total variability of artifact-adjusted RR heartbeat interval series within a specified time interval. The standard deviation of the RR distances is calculated in milliseconds (ms). In an exemplary analysis sequence according to Figs. 19-21, a comparison of concurrently detected HFI, heart rate, and SDNNRR events is performed. The event theory of the SDNNRR curve according to FIG. 21 associated with the preceding FIGS. 19 and 20 in detail:
19, a multiplicity of primary significances 1-n has again been detected, which with substantially simultaneously occurring secondary significances 1'-n 'in the heart rate progression diagram according to FIG. 20 as well as with respectively corresponding secondary significances T'-n "in FIG SDNNRR waveform diagram according to Fig.21 coincide in time. This is followed by a one-time comparison based on the previously performed graphical evaluation regarding extrasystoles and QPA. These are transferred to the SDNNRR data set and compared within the time intervals T1-T11.
The time interval T1 shows (almost) no LFD and thus a negligible HFl-induced influence.
The SDNIMrr is increased due to the increased variability from the heart rate. During undisturbed sleep, history fluctuation within the natural range, accompanied by a QPA in the optimal 4: 1 ratio, shows a high RSA, which shows in the form of the logRSA progression (second curve below the SDNNRR curve and a flat, lowered heart rate In practice, the SDNNrr and logRSA traces and the trapped areas between corresponding trajectories and abscissa corresponding to the time axis are color coded, eg, blue and violet, but this color coding was again resolved into a black and white dot graph in Fig. 21 Slash-hatched areas (= high RSA) of the subsurface logRSA curve are associated with deep sleep phases, with a deep sleep phase just ending in time interval T1, while a new deep sleep phase begins between time intervals T2 and T3 (again hatched), but by the in Fig.19 ersic can no longer properly form HFI load or LFD tips.
Between time intervals T2 and T6, typical staircase courses in the LFD (shown in Fig.19 as primary significances 1a, 2a, 3a) with significant peaks at the beginning of the staircase form and then temporarily (approximately 15 minutes) of constant LFD are shown. These LFD peaks (shown in Fig.19 as primary significances 1, 2, 3, 4) result in coincidental extrasystoles in the heart rate progression according to Fig.20 (plotted as secondary significances T, 2 ', 3', 4 '), respectively with increased autonomous activation in the SDNNRR. At the same time, in the phases of constant LFD, the lowering of the QPA and a corresponding reduction in the SDNNRR (shown in FIG. 21 as secondary significances 1a ", 2a", 3a "). Also noteworthy are analogous flattening of the heart rate curve (shown in FIG. 20 as secondary significances 1a ', 2a', 3a ') in each case shortly after the end of the extrasystoles 1, 2, 3 marked in FIG. 19 or, respectively, during the horizontal segments of the denoted LFD. Stair shape (= primary significances 1a, 2a, 3a).
The time intervals T2-T6 have been marked with a horizontal black line as the 1st (activation) level in the SDNNRR progress diagram according to FIG. In time interval T6, the duration of the constant LFD (FIG. 19) is shortened, which is accompanied by an increase in the SDNNRR (increased activation) and decreases again to the first level by a longer lowered LFD in time interval T7. In time interval T8, there is an accumulation of primary significances (marked by arrows in FIG. 19) in the form of extreme short changes, i. Increase and decrease of the LFD with many significant LFD peaks, which according to Figure 20 result in corresponding extrasystoles and increase the QPA significantly. In the SDNNRR course, secondary significance is shown by the strong increase in autonomic activation (see indicator arrows in Fig. 21 and an oblique line drawn between the 1st and 2nd levels, respectively, to illustrate the increase in SDNNRR).
In time interval T9, the LFD drops again to a very low level according to FIG. 19 and a deep sleep phase begins again (see hatched area in the logRSA curve). At the same time, in this time interval with the onset of the deep sleep phase, a compensation initiated by the autonomic organ system for the preceding activation sets in, which shows itself in the heart rate progression by numerous interval extensions (vertical thin lines below the thick average curve). In the SDNNRR, this initial activation is shown with a lowered heart rate, optimal QPA 4: 1, and high RSA (hatched area). The SDNNRR reduces until the next increase of the LFD (FIG. 19) to a level 2 plotted in FIG. 21 and thus differs completely from an externally caused activation. With time interval T10 then begins a massive increase and chaotic course of HFI (Fig.19) with numerous LFD peaks, which exerts the strongest stimulus on the organism and the SDNNRR increases simultaneously to the highest (3rd) level. Even a short-term decrease of the LFD according to Fig. 19 does not lower the autonomic activation of the organism below the third level.
Especially in the present time interval T10, it becomes clear how the test subject's organism tries to maintain the deep sleep phase, visible in the taut area in the logRSA course and the heart rate with numerous vertical lines, but this internal compensation attempt becomes due to the primary significance of a rapidly changing LFD with numerous needle-shaped peaks severely affected or chaotized. The organism considered in the previous analysis sequence shows phases such as under day activation and a dissolving sleep architecture.
In particular, a heart rate variability (HRV) frequency information data set that can be represented as a "spectrogram" can be analyzed as the HRV-associated autonomic functional parameter. The spectrogram already explained above represents activation over a frequency bandwidth of substantially 0 to 0.5 Hz, each associated with a specific frequency range of vegetative organ functions of the test subject. An analysis sequence as illustrated in FIGS. 9-11 will now be described. In this case, a comparison of the HFI curve according to FIG. 9, of the spectrogram according to FIG. 10 and of the VQ curve according to FIG. 11, which has been calibrated with respect to the time intervals T1-T11, is undertaken.
The loss of the sleep structure of the test subject already mentioned in the preceding evaluation of the SDNNRR course, according to high-frequency immissions, can hereby be represented significantly. It should be noted in advance that the interaction of breath, blood pressure and blood flow is an expression of the vegetative control of the heartbeat (VQ: LF / HF), whereby the vagus (HF: rest and recovery) in particular causes the modulation of the respiration (RSA) and Blood pressure and circulation of sympathetic control (tension) subject.
According to Figure 9 shows in the time interval T1 a very low and therefore negligible LFD.Entsprechend a good sleep architecture shows a deep sleep phase (plotted arcuate line / indication P1) with a pronounced RSA (frequency range between 0.2 -0.3Hz). In the BRAC (basic rest-activity cycle), which is also already defined at the beginning, the sinusoidal or convex course of the RSA (= deep sleep phase) alternates with the REM phase. This change usually takes place with 90-120 minutes deep sleep phase (recovery phase) and 5-20 minutes REM phase (activity phase, the sinusoidal course of the BRAC is also shown in Fig. 27 in an exemplary manner). In the VQ, there is already a clear reduction of the activity level with a still slightly dominant sympathetic nervous system (ratio LF / HF = 2: 1).
The color codings blue, white, yellow, orange, red used in the following exemplary embodiment (each frequency-specific pixel shown in the spectrogram also has depth information relating to activation or amplitude strength - see also the initially defined assignment) according to the three-dimensional information content of the spectrogram correspond to one of the Practical original evaluation. The spectrograms of FIGS. 10, 22 and 26, which are converted into monochrome or pure black-and-white pixeling, can give only a rough impression of this color gradation.
The event history relating to FIG. 10: In the time interval T1 of the spectrogram, predominantly blue (in FIG. 10: dark) with some white amplitudes show up, which still points to an irregular inclination. In contrast to the REM phase, there are no red amplitudes (in the frequency range between = and 0.1 Hz) concerning peripheral tissue perfusion (0.017 Hz) or activation of the neuronal system with red amplitudes until just below 0.1 Hz. There is also no strong activation in the Range around 0.1 Hz, which corresponds to the blood pressure range, but usually also in the REM phase is omitted and only in the daily rhythm (waking phase).
Time interval T2 lies between the deep sleep phases P1 and P2, the REMPhase drawn in FIG. 10 coinciding in the spectrogram with a first increase in the LFD (see HFI profile according to FIG. As a result, the circulation starts and in the VQ there is an increase in the activity level (sympathetic), which is the normal course of the change between deep sleep and REM phases. However, the HFI trace shows a sharp change, i. an increase and decrease of the LFD and some acicular tips (primary significances, not indicated here by Pos. No.) which are detectable as an irritant to the autonomic system. This is also evident in an unnaturally strong increase in the VQ (above 10: 1) in FIG. 11 and in the blood pressure activation (red amplitude) just above 0.1 Hz in the spectrogram according to FIG. Such a blood pressure activation would normally not arise even in the REM phase and is therefore already detected as a secondary significance in the method according to the invention. This becomes clear in comparison to the second REM phase in the time interval T6, in which the HFI curve again shows a staircase form with a constant (approximately 10 min.) Constant LFD, which exerts less of a stimulus on the organism. There is no blood pressure activation in the frequency range of about 0.1 Hz, the VQ does not rise so high (VQ = 5: 1), the general activation in the spectrogram according to FIG. 10 also falls in the frequency range below 0.1 Hubstantially lower.
Also worth mentioning are three small peaks occurring in the time interval T6 in the VQ progression diagram, which reach an LF / HF level of approximately 10: 1. These peaks or needle-shaped elevations are detected as secondary significances (indicated in Figure 11 by three black dots without Pos.No.). As shown by a comparison with the event history in the HFI progress diagram, these secondary peaks coincide in time with primary significances (acicular tips or striking endpoints of LFD staircase forms) also marked by three dots in FIG. 9 and are therefore evaluated as HFI-induced. This example shows how even relatively small LFD peaks or primary significances can possibly lead to significant reactions in HRV-FP curves.
The time intervals T2 to T5 according to FIG. 10 show the course of a second deep sleep phase P2. The RSA is clearly formed in the spectrogram at the beginning, evident from the white amplitudes at the beginning of the sinusoidal or arcuate course (in the black-and-white resolution according to FIG. 10 only recognizable), but then decreases successively in time coincidence with the entrance of several Although LFD staircase progression of the LFD hereof activates the organism (sympathetic) less than described above, the regular time intervals of the LFD peaks or primary significances, for example, act like LFD peaks or primary significances 1, 2, 3, 4 according to FIG a meta-rhythm that attenuates the respiratory rhythm in heart rate (RSA). This becomes clear again in contrast to time interval T10, in which the chaotic HFF course with numerous LFD tips activates the organism more strongly, but the RSA, although fragmented, has remained in its intensity as it did at the beginning of the deep sleep phase P4.
This also corresponds to the above-described evaluation of the SDNNRR in the KurvelogRSA (FIG. 21): The shaded regions of the logRSA path coincide with the pronounced RSA, while the logRSA clearly decreases during the periods of attenuated RSA. The consequence of the attenuation of the RSA is then seen at the beginning of the deep sleep phase P3 in the time interval T7 (FIG. 10), in which the RSA is barely recognizable and despite low LFD in the HFI curve according to FIG. 9, only in the VQ (FIG shows a below-average reduction, especially since the tone of the vagus is weakened.
The VQ in the deep sleep phase P2 shows a largely lowered course over the time intervals T2-T5. However, VQ depression is massively disrupted four times with spontaneous activation peaks to the 10: 1 ratio (plotted in Figure 11 as secondary significances 1 ", 2", 3 ", and 4"), which coincide completely with LFD peaks in HFI progression 9 (shown as primary significances 1,2, 3 and 4).
At the same time, very strong activation at the onset of blood flow up to blood pressure activation (seen as red-yellow-white color coding in the spectrogram; in the present black and white pixel graphic of Figure 10, only as white columns in P2) The areas of these - and also the other - bright pixel clusters which are close to the abscissa or between 0 and 0.05 Hz are respectively marked red in the original spectrogram and color corresponding to the highest activation, but this red coding occurs in the present FIG. As a result of the black-and-white transformation of the graphic, it was also converted into a predominantly black pixel accumulation like the oppositional blue regions, which indicate a minimal or nonexistent activation of respective organ systems).
As can be seen in the VQ progression according to FIG. 11, the vagus nerve, which has tried to compensate, can still interspersed again and again between the aforementioned LFD peaks or primary significances 1-4 (step-shaped course in the VQ diagram / FIG. 11), as shown in FIG Fig. 10 even results in a deepening / deepening improvement compared to the previous deep sleep phase P1 with significantly less white coded amplitudes. Due to the four-time significant interruptions in the time intervals T3-T5, however, the organism does not reach such a deep depression of the sympathetic that the vagus would predominate and correspond to an intact sleep architecture.
A third deep sleep phase P3 beginning in the time interval T6 (after the second REM phase) exhibits almost no RSA activation in the spectrogram according to FIG. 10 (recognizable by the dark range between 0.2 and 0.3 Hz). As described above, the RSA has been attenuated by the HFi-induced, staircase-shaped metarhythm, and after the REM phase can not be re-amplified in similar LFD staircase patterns or primary significances in the HFI course. Although in the time interval T7 a lowering of the LFD to be recorded in the HFI curve is associated with a decrease in the VQ in FIG. 11, this subsidence remains below average. In this deep sleep phase P3, the time course, which would have continued due to the BRAC cycle explained at the outset (a deep sleep phase usually lasts 90-120 min.), Abruptly by the very strong increase of the LFD beginning in the time interval T8 or further in FIG apparent needle-shaped significances were aborted. In the spectrogram according to FIG. 10, activation of the frequency range assigned to the perfusion (0-0.1 Hz) and, in the previous event history, the strongest activation of the frequency range assigned to the blood pressure (range around 0.1 Hz) are once again shown.
Shortly after 4 o'clock (time interval T8), a short-term low-LFD is again shown in the HFI curve according to FIG. 9, which is accompanied by a simultaneous decrease in the VQ in FIG. 11 and the transition to a fourth deep sleep phase P4. At this point it should be noted that not only elevations, but also pre-defined subsidence, thus temporary "improvements" of the autonomic balances as secondary significances in the HRV-FP course, e.g. can be detected in the VQ curve if a correlation is given with detected or characteristically defined HFI events / start-up sections. The exact location of the following three REM phases is now only more difficult to recognize from the spectrogram since, at the same time as the transition between deep sleep phase P3 and P4 (change from time interval T8 to T9), a strong LFD peak or primary significance in the HFI curve according to FIG Fig.9auftritt. In the time interval T9, a longer time span (almost) without LFD is shown for the first time, thus simultaneously the development of a more pronounced RSA in the spectrogram (FIG. 10) and a significant decrease of the VQ until the first prevalence of the vagus in FIG. In the fourth low-pass phase P4, there is again a normal course of the sleep architecture, without activation of the frequency ranges blood circulation or blood pressure, with predominantly blue coded, i.e., weak amplitudes. However, the deep sleep phase P4 abruptly disturbed from the time interval T10.Beginnt in the spectrogram at about 04.30 clock with significantly pronounced RSA without red, but predominantly only blue coded amplitudes (in terms of color coding or their graphic representation still applies the foregoing), so has the Deep sleep phase only 45 min. from the usual 90-120 min. (BRAC) when they suddenly show a strong increase in LFD in the HFI course at around 5:15 pm (seen as a vertical, comb-shaped primary significance in Figure 9), as a result of strong activation of the sympathetic massive influx of blood for activation of the blood pressure (visible in the columnar activation in the spectrogram detected as secondary significance, the activation in this case exceeds over 70% of the entire frequency bandwidth of the spectrogram or in the purely mathematically evaluable H RV-F frequency-1, as shown in FIG nformationsdatensatz).
This strong activation of the entire organism then stops at the same time as the chaotic sequence of LFD peaks or primary significances in the HFI curve according to FIG. 9 to 07.15 o'clock. In the spectrogram, the RSA / frequency range 0.2-0.3 Hz) in their previous arcuate course is completely ruptured, the course of the BRAC (change of deep sleep and REM phases) is completely lost. Only at 07.15 clock in time interval T11, as it was for about 20 min. According to the spectrogram, a fifth deep sleep phase P5 can be established - with short exposure to blood flow, no blood pressure activation (indicated by the blue or dark coding of the pixels located in the respective frequency ranges) and one shown in FIG Lowering the VQ with prevalence of the vagus.
For the evaluation of the HRV frequency information data set or spectrogram or for recognition as secondary significance, suitable criteria or threshold values can be defined. For example, it can be evaluated as secondary significance if the amplitude strength of respective frequency ranges exceeds 20%, 30%, 50%, 70% or 100% increased.
In an analysis variant shown with reference to FIGS. 22-29, the RSA which can be represented in the HRV frequency information data set or spectrogram or its course lying in the frequency range from 0.2 to 0.3 Hz is considered, with a substantially convex toward the upper end of the frequency bandwidth vaulted course of activated organ functions, indicated in Figure 22 by white arcs 1'-5 ', is evaluated as HFI-induced secondary significance. The convex curves or secondary significances 1'-5 'extend over time intervals between 30 and 120 minutes and also show a decreasing activation or amplitude strength (only insufficiently recognizable in the present FIG. 22 due to black-and-white transformation, but in the original spectrogram In the case of significant HFI loads, the previously described chronobiological rhythm of the BRAC (basic rest-activity cycle) is often reflected in the subject RSA frequency range which indicates activation by the sympathetic and parasympathetic - and thus not by the respiratory center. In sleep, the BRAC manifests itself as a change between deep sleep and REM phase (see also black, arcuate indication lines in Fig. 24 and Fig. 25), with the REM phases also showing on the lowered arch forms and end portions of the convex RSA progress sections, respectively ,
An apparent in Fig.22 mirroring the BRAC in the otherwise essentially barkenförmi¬gen or parallel to the time axis - in a frequency range of usually about 0.25 Hz - running RSA course during sleep (see Fig.26) is therefore considered a suitable criterion for detecting an HFI-induced disturbance of the vegetative regulatory mechanisms of the test subject. An increase in RSA progression associated with weakening of the respiratory center, e.g. at approximately 0.3 Hz, an anti-cyclic slowing of the respiration indicates a decrease in RSA progression indicates anticyclic acceleration of the test subject's breath and thus an increase in sympathetic activation; in total, the respiratory rhythm is at least partially chaotized. Corresponding anomalies or secondary significances are also detectable in the simultaneously considered QPA course (see FIG. 25).
The present analysis of the spectrogram is also an example of the possibility of a time-delayed occurrence or detection of secondary significances, especially since the illustrated HF1-related RSA / BRAC profile can also be observed in a period in which the test subject does not have any relevant HFI. Is no longer exposed to stress or primary significance. FIG. 22 already shows such a scenario - a sleep phase of the test subject following approx. 24 hours of strong HFI exposure was shown, whereby no relevant LFD was measured on this second night of the IT dosimeter. Only later, in this case in the third HFI exposure-free night shown in FIGS. 26-27, does the described pathological tendency in the RSA frequency range disappear and the RSA curve smoothes again or shows an approximately bar-shaped form running parallel to the time axis (see Fig.26).
Since the above-described arrhythmias or secondary significances in one or more HRV-FPs are accompanied by relevant physiological impairments or a chaoticisation of the vegetative balance, upon detection of a defined number of such significances per reference time interval, a positive assessment of the presence of a human biologically relevant influencing of the vegetative regulatory mechanisms of the test subject. With prolonged nocturnal disturbances, then, in addition to reduced vagal interference, even during the day, loss of cardiac vagal modulation of the heart rate with reduced performance occurs. Sleep-related respiratory disorders (e.g., hypopneas) are also associated with symptoms such as falling asleep and staying asleep, restlessness, nightmares, depression, and anxiety. As a further consequence of a Meiatonindefizit according to predominantly Sympathikusak¬tivierung is recorded, a disproportionate Kortisolbildung in specific brain areas and a sinking memory in declarative memory contents.
FIGS. 30-34 illustrate a particularly advantageous variant for evaluating the HRV-FP and HFI courses. Fig. 30 exemplarily shows an HFI / LFD waveform, Fig. 32 shows a heart rate trace detected at the same reference time interval. Instead of the heart rate trace, any other HRV-FP trace may also be analyzed in a manner described below.
As previously described, the HFI and HRV-FP measurement / computation values obtained from two different signal channels were scaled (i.e., maximum value selected) at the same number of sub-time intervals in each of the 7-second sub-time intervals. First, both the measurement or calculation values constituting the HFI waveform and the HRV-FP trace are converted to a chronological gradient sequence GHrv-fp and GHfium, respectively, which are composed of a plurality of successive gradients, i. Difference formations between two successive HFI or HRV-FP measurement or calculation values: gradient = signal (index) - signal (index + 1). Fig. 31 shows how absolute values of the gradients are displayed in the form of bar graphs. The analysis of the gradient sequences GHrv-fp and GHfi as described below with regard to their correlation can be made on the basis of the graphical processes or also in a purely computational manner. The gradients of the gradient sequences GHrv-fp, Ghfi can also be defined and evaluated as vector variables.
The first gradient sequence Ghrv-fp corresponding to the HRV-FP profile, in the present case the heart rate, and the second gradient sequence GHfi corresponding to the HFI profile are arranged on both sides of a mirror axis extending parallel to the time axis or corresponding thereto. The individual, bar-shaped, or gradient-shaped gradients of the gradient sequences GHRv-Fp and GH respectively depend on the mirror axis labeled zero (in this case horizontal) and are orthogonal to it. In the present exemplary embodiment according to FIG. 31, the gradients of the first gradient sequence Ghrv-fp rise more or less upward according to their absolute value, while the gradients of the second gradient sequence GHfi extend more or less far downwards according to their absolute value. Regarding the graphic arrangement and evaluation of the gradient sequences GHrv-fp > GHfi are once again ready for a multitude of possibilities. In principle, all the methods of visualization and analysis already mentioned above for the analysis of the HFI and HRV-FP profiles can also be used here. Thus, some of the starting point of the gradient sequences GHrv-fp, Ghfi or their zero level may also be spaced apart from the mirror axis, the diagrams could run in the vertical direction, be color coded etc.
For better visual comparability, the gradient sequences GHrv-fp and GHfi may be increased or decreased in terms of their extension normal to the time axis. In the example according to FIG. 31 (also in FIGS. 33 and 34), the size of the gradient sequence GHr corresponding to the HFI curve or the extent of the HFI gradients normal to the time axis, thus here in the vertical direction, has been increased. Such an approximation or scaling is preferably carried out with reference to average values of the HRV-FP and HFI courses or their gradient sequences Ghrv-fp and GHfi within the respectively considered reference time interval.
Fig.33 shows how to the graphs of the gradient sequences Ghrv-fp. GHfi envelopes (trend curves) were created. In this way, any correlation in the time course of the gradient sequences GHrv-fp. GHfi or the simultaneous occurrence of short-term primary and secondary significances can be shown particularly clearly.
In an alternative representation variant according to FIG. 34, the gradient sequences GHrv-fp and GHfi or their envelopes are shown overlapping each other so that any correlations between the two gradient sequences GHrv-fp and Ghfi or HRV-FP and HFI progressions are made directly visible. As can be seen in particular in FIG. 34 (see also FIGS. 31 and 33), the gradient sequences GHRv-Fp and GHFi run over long paths approximately synchronously or at least with a conforming tendency. Significant maxima in the first gradient sequence GHrv-fp coincide, in terms of their time of occurrence, with corresponding maxima in the second gradient sequence GHfi representing the RFI responses to the test subject.
The envelopes shown in FIG. 33 and FIG. 34 each affect the tips of the gradients facing away from the zero level or the time axis, respectively, or are nestled against them. Herein, for example, virtual tangents to which the envelope approximates asymptotically or finds a reversal point may be applied to the designated gradient peaks. In the present evaluation, the gradients of the gradient sequences GHRv-Fp and GHfi were computationally generated by Nadaraya-Watson kernel regression.
Between the gradient sequence GHrv-fp corresponding to the HRV-FP profile and the gradient sequence (GHFi) corresponding to the HFI or LFD curve, a correlation value, e.g. is determined as a weighted cross-correlation coefficient. The correlation value may e.g. be specified as a decimal value or even in percent. For example, in the present evaluation, it has been determined that two time-correlating HRV-FP and HFI gradients, i. primary and secondary significances have a correlation of 1, while uncorrelated HRV-FP and HFI gradients are assigned a correlation of 0, and exactly opposite HRV-FP and HFI gradients correlate -1. In the present evaluation according to FIGS. 30-34, a correlation of 0.656 was achieved in this way, which is to be regarded as a noteworthy example for the preponderance of a human biologically relevant influencing of the vegetative regulatory mechanisms of the test subject.
In a special evaluation on QPA / RSA / respiratory rate shown in FIGS. 16-18, the physiological role of these functional parameters or the HF1-related influence on sleep is illustrated. As already mentioned, the oscillation of the RSA in the spectrogram is found in a healthy adult in the frequency range between 0.2 and 0.3 Hz. With the onset of RSA, a respiratory-synchronous change in the heart rate is shown in the heart rate, which makes the RR intervals shorter when inhaled ( accelerated) and prolonged (exhaled) during exhalation. It is mainly mediated by the parasympathetic / vagus, which becomes active when exhaling. In addition to the spectrogram, RSA is also shown in logRSA, which accurately calculates this rapid respiratory-induced change in RR interval changes and thus reflects the rate of cardiac rhythm modulation by respiration (compare also the hatched areas of the logRSA of Fig. 21 with the RSA in the spectrogram) Figure 10).
The function of RSA is to optimize the gas exchange on each respiratory cycle by timing the pulmonary circulation heartbeats to the airway ventilation (through the breaths). The strength of the RSA depends on: tone of the parasympathetic nervous system (modulation capacity of the nervous system), respiratory rate, tidal volume.
Has the depression of the tidal volume only a minor effect on the amplitude or
Strength of RSA and the direct influence of radio frequency on the nervous system is not representable here, so a change in the respiratory rate directly affects the RSA and is also in a certain dependence on functional requirements (such as hypopneas). An increase in the rate of respiration of e.g. 11 to 15 breaths per minute increases the frequency of the RSA in the spectrogram from 0.19 Hz to 0.25 Hz and vice versa, thus showing the influence of the respiratory rate on the heart rate. A maximum of the RSA amplitude should be observed between 12-15 breaths per minute. Unlike RSA, which allows heartbeats to be beat in a particular rhythm during a respiratory cycle with a vagus, in QPA the number of heart beats optimizes to one (4: 1) respiratory cycle, thereby maximizing the vagal impact of RSA. This relationship is shown graphically by the amplitude in the spectrogram and mathematically following model calculation. Thus, RSA and QPA together optimize the performance of the heart and lungs during sleep, with strong coupling (4: 1) with a shortest possible RR interval during a late inspiration (which accelerates the heart rate) and a longitudinal RR interval immediately before the beginning of the next Inspiration (has slowed the heart rate again) shows. A wave motion is formed which supports maximum acceleration of the heart rate in the inspiration. This vegetative synchronization relieves the heart and thus provides energy for the physical regeneration. The effects of the primary significances or HFI / LFD peaks have so far been clearly evident in the temporal comparison with the triggering of extrasystoles and the sudden increase in QPA progression (see FIGS. 12 and 13 and FIGS. 16 and 18). This sudden increase in QPA by external activation thus destroys the ideal coupling (= 4: 1) between heartbeat and respiration occurring during sleep, depending on the HFI progression or temporal incidence of the primary significances.
In the QPA course according to FIG. 18 (see also time interval scale according to FIG. 15), the optimal integer ratio of 4: 1 is shown in the first time interval T1, and the organism always shows it with every longer decrease (time intervals T4, T7, T9) of the LFD and that it increases with each LFD and increases the respiratory rate. This episode has not yet been discussed and is reflected in the spectrogram as a result of the effect of altered RSA in the HRV analysis. This mathematical and graphical relationship is explained below: The mathematical calculation is based on an HRV diagram sequence shown in FIGS. 16-19, from which, for example, nine comparison points were selected over all time intervals which always have a primary significance in the form of an LFD peak (FIG. 2, 3, 4, 6, 7, 8, 9) or coincide with an almost LFD-free HFI progression phase (1, 5) and in each case with the simultaneously determined heart rate according to FIG. 17 and the QPA according to FIG were compared.
From the above calculation, it is also possible to derive a simple high-level calculation for the additional load of the heart, which is exactly verifiable from the original measurement data of the ECG: The average number of heartbeats of the comparison points 1'-9 'marked in FIG. 17 is 67.33 / minute In comparison, the heart rate in the first, largely LFD-free LFD time interval T1 (see comparison point / primary significance 1 in FIG. 16) is 53 / min. This results in an additional burden of 6,878 strokes at eight hours of sleep. With intact sleep architecture, the heart rate only slightly increases in the transition between deep sleep phases and REM phases, which is why this was neglected in the calculation.
The corresponding spectrogram according to FIG. 10 shows clearly pronounced white amplitudes of the RSA with sinusoidal or convex profile and strong coupling in the time intervals T1 and T9 (corresponding to the calculated comparison points 1 and 5 according to FIG. 16). Time interval T10, on the other hand, clearly shows, at the white, chaotically structured amplitudes (indicated by a black wavy line in Fig. 10), a fragmentation of the RSA in the range between 0.2 - 0.3 Hz (RSA / BRAC) and beyond, manifesting itself in represents the above calculation. Here is the correspondence between respiratory rate and RSA. It practically indicates the above calculated respiratory rate. The objective relationship between primary significances in the HFI course and the simultaneous change in the respiratory rate or its chaotic course, which acts via the RSA on the respiratory synchronous heart rate, can therefore be used in a preferred analysis variant of the method according to the invention as a criterion for the detection of HF1-induced secondary significances.
In general, any external perturbation intermittent to physiological pulse-breathing coupling proves to be. any disproportionate acceleration of the heartbeat in relation to the respiratory rhythm as a stress factor, which does not necessarily have to assume pathological proportions, but is at least perceived as subliminal dyspnea (respiratory distress). A, albeit subconscious, respiratory distress or oxygen deficiency always causes a subconscious, anxiety reaction in humans on a vegetative level. The increased incidence of anxiety disorders associated with fatigue / depression symptoms, as seen in clinical psychology, is certainly recognized as a multicausal phenomenon, but may be given diagnostic and preventive support by the HRV / HFI analysis method of the present invention.

权利要求:
Claims (30)
[1]
A method for analyzing the interaction of electromagnetic high-frequency immissions (HFt) with vegetative regulatory mechanisms of a test subject, comprising the following steps: Determining the course of an intensity of radio-frequency immissions (HFI), preferably measured as power flux density, with respect to the test subject - the determination of the variance in the time interval of successive heartbeats recorded by means of electrocardiogram measurement (ECG) and the progression of their determined, associated with the heart rate variability (HRV) vegetative function parameters (HRV-FP) of the test subject over a defined period of time; HFI and HRV-FP profiles are at least partially synchronized and at least one common reference time interval (T | _n) calibrated, - performing an analysis of the HFI curve, in which primary significances (1-n), ie temporary, i. within the reference time interval OY ,,) detects time-limited deviations of the HFI curve from defined static or dynamic mean values or HFI reference data, in particular significant increases and / or decreases in the HFI curve with respect to temporally preceding and / or subsequent sections of the HFI curve - Perform an analysis of the HRV-FP history in which secondary significances (Tn '), temporary, ie within the reference time interval OYn), time-limited deviations of the HRV-FP curve from HRV-FP reference data or rule traces derived from empirical data inventory and / or generated dynamically, i. derived from current measurement / calculation data or HRV-FP progressions, in particular significant increases and / or decreases in the HRV-FP progression with respect to time-preceding and / or subsequent sections of the HRV-FP run; upon detection of a defined number of coincident, i. in the HFI course and in the HRV-FP course, a positive evaluation of the presence of a human biologically relevant influencing of the vegetative regulatory mechanisms of the test subject takes place in a substantially coincidental primary and secondary significance (1-n, 1'-n ').
[2]
2. A method for analyzing the interaction of electromagnetic high-frequency immissions (HFI) with vegetative regulatory mechanisms of a test subject, comprising the following procedural steps: Determining the variance of the time interval of by means of electrocardiogram measuring method (ECG) detected, successive heartbeats and the course thereof determined, with the heart rate variability (HRV) associated vegetative function parameters (HRV-FP) of the test subject over a defined period, - comparing the currently determined HRV-FP course with characteristic, HFl-induced secondary significances or with HRV-FP reference data or Regulating processes, which respectively ausaktueller experimental design or from at least one of these preliminary Versuchsanordnnung with the current test subject or third test subjects were determined and now in the form of digital or analog data, preferably in the form of graphic gradients, alphanumerischen or algorithmic data, functions or function graphs are kept on a storage device; - Upon detection of a defined, at least partially present agreement of the currently determined HRV-FP curve with characteristic, HFl-induced secondary significances or deviation of the currently determined HRV-FP curve from the HRV-FP reference data, a positive assessment of the presence of a human biologically relevant Influencing the vegetative regulatory mechanisms of the test subject takes place.
[3]
3. The method according to any one of claims 1 or 2, characterized in that as a vegetative function parameter (HRV-FP) as a "spectrogram " variable heart rate variability (HRV) frequency information record is analyzed which reflects the activation across a defined frequency bandwidth of substantially 0 to 0.5 Hz each to a specific frequency range of associated vegetative organ functions of the test subject, such as: 0.04 to 0.15 Hz: Low Frequency (LF) -with correspondence: predominantly sympathetic activity, to a lesser extent also vagal activity, assignment in particular of the blood pressure and circulation rhythm, 0.15 to 0.40 Hz: High Frequency (HF) with correspondence: vagus activity; Assignment of, in particular, respiratory functions, preferably the respiratory sinus arrhythmia (RSA) reflecting the modulation of the cardiac rhythm by respiration, and activation of respective frequency ranges or autonomic organ functions can be taken from the HRV frequency information data set in the form of amplitude intensities preferably visualized by color coding is detected as HFl-induced secondary significance (Tn ') when spontaneous activation of previously unavailable or negligible intensity activated or autonomic organ functions over a bandwidth of at least 0.05 Hz, preferably at least 0.1 Hz , Particularly preferably over a Band¬ width of more than 0.2 Hz is carried out.
[4]
4. The method according to claim 3, characterized in that it is detected as a secondary significance (1'-n '), if a short-term, preferably not more than ten minutes, more preferably not lasting more than one minute activation of before still not or invernachlässigbarer intensity of activated frequency ranges or vegetative organ functions takes place, which is in the case of a graphical visualization of the HRV frequency information data set as approximately needle-shaped, orthogonal to the time axis survey, in a preferred method variant as a sufficient criterion for detection as se Secondary significance (1 '-n') must be a simultaneous activation of more than 50%, preferably more than 70% of the total, substantially 0 to 0.5 Hz frequency bandwidth of the HRV frequency information data set or spectrogram.
[5]
5. The method according to claim 3 or 4, characterized in that it is evaluated as activation or secondary significance (1'-n '), if the amplitude strength of respective Frequenzbe¬reiche or vegetative organ functions by more than 20%, preferably by is increased more than 50%, and the analysis method may also include the possibility of weighting detected secondary significances (1'-n ') on the basis of respectively determined amplitude intensities, which are color-coded in the case of graphical visualization, and / or on the basis of the size of the respectively activated frequency bandwidth , where events or secondary significances (Tn ') at which high amplitude strengths or activations reaching over broad ranges of the frequency bandwidth of the HRV frequency information data set have been determined, with a larger factor for the assessment of the presence of human biological relevant influence on the vegetative test subject regulation mec be charged as those events. Secondary Significancies (T-n '), in which lower amplitude strengths or over a smaller areas of the frequency bandwidth of the HRV frequency information record reaching activations were determined.
[6]
6. The method according to any one of claims 1 or 2, characterized in that is analyzed as a vegetative function parameter (HRV-FP) representable as a "spectrogram" HRV frequency information data set, which over a defined frequency bandwidth of substantially 0 to 0.5 Hz reproducing a specific frequency range of associated autonomic organ functions of the test subject, wherein a frequency range of essentially 0.2 to 0.3 Hz of the HRV frequency information data set lying temporal course of vegetative organ functions is analyzed, which reflects the respiration-induced modulation of the heart rhythm respiratory sinus arrhythmia (RSA), wherein a convexly curved course towards the upper end of the frequency bandwidth is evaluated as HF1-induced secondary significance (1'-n '), the convexly curved course preferably over a time interval between 30 and 160 M. In particular, preferably between 60 and 120 minutes extends.
[7]
7. The method according to any one of claims 1 or 2, characterized in that as vegetative function parameters (HRV-FP) the time course of the heart rate of the test subject is analyzed, and thereby the occurrence of short-term arrhythmias over the temporally preceding and following heart rate course as a secondary significance (1 '-n') is detected.
[8]
8. Method according to claim 7, characterized in that the arrhythmia of the heart rate and thus as secondary significance (1'-n ') the occurrence of extrasystoles is detected, whereby an extrasystole in the course of the in the unit: [number of heartbeats per unit of time, vorzugs¬weise per minute] measured heart rate as a short-term, erratic increase in the heart rate over the respective chronologically preceding and subsequent heart rate course and thus as a significant shortening of the interval between two consecutive heartbeat, in a preferred process variant, erratic increases in the heart rate can then be recognized as extrasystoles, if these increases are in each case more than 30%, preferably more than 50%, particularly preferably more than 100% of the chronologically immediately preceding or following heart rate or in relation to the average heart rate determined in the respective measurement period.
[9]
9. The method according to any one of claims 1 or 2, characterized in that the vegetative function parameter (HRV-FP) is analyzed as the time course of respiratory sinus rhythmic respiratory sinus arrhythmia (RSA), which is preferably converted to logRSA by means of a decadic logarithm.
[10]
10. The method according to any one of claims 1 or 2, characterized in that as a vegetative function parameter (HRV-FP) the time course of the pulse-breath quotient (QPA) of the test subject is analyzed, preferably those events as secondary significances (Tn ' ) are detected, in which the QPA curve deviates from a ratio pulse / breath = 4: 1 by more than 20% and / or in which temporary peaks or maximum values of the QPA profile occur.
[11]
11. The method according to any one of claims 1 or 2, characterized in that as a vegetative function parameter (HRV-FP) the time course of the vegetative quotient (VQ) of the test subject is analyzed, which is the ratio of predominantly the sympathetic, in particular the Blood pressure rhythm associated low-frequency, substantially between 0.04 - 0.15 Hz extending HRV frequency ranges (LF) to the parasympathetic, insbesonde¬re the respiratory modulation associated high-frequency, substantially between 0.15 and 0.40Hz extending HRV frequency ranges (HF) results, where as secondary Significances (1'-n ') are detected: temporary increase of the vegetative quotient (VQ) by more than 20% (compared to temporally immediately preceding VQ values or to a VQ mean value determined in the reference time interval) and / or transgressions or Sympathetic dominance of the vegetative quotient of more than 3: 1, preferably more than 5: 1.
[12]
12. The method according to any one of claims 1 or 2, characterized in that as a vegetative function parameter (HRV-FP), the time course of the SDNNRR (standard deviation of normai-to-normal intervals) or a SDNN corresponding statistical characteristic such. SDNNIDX / ASDNN of the test subject is analyzed, i.e. an HRV-associated statistical spread around the mean of the heartbeat interval duration or its differences representing the change in total variability of artifact-adjusted RR heartbeat interval series within a predetermined time interval comprising substantially all frequency ranges.
[13]
13. Method according to one of claims 1 to 12, characterized in that as secondary significances (1'-n ') in one of the above-mentioned HRV-FP curves, all events or run segments are detected which, when graphically represented, have a slope k > 1 or a pitch angle α of more than 45 °, preferably of more than 70 °, more preferably of approximately 90 ° with respect to the time axis (slope k- »°°) or an approximately needle-shaped elevation of the HRV-FP course the above-mentioned events or detection criteria must each exceed at least 20%, preferably at least 30%, of the average HRV-FP profile.
[14]
14. The method according to any one of claims 1 to 13, characterized in that from the time points of a detection of HFl-induced secondary significances (1'-n ') each time at least one sub-time interval is opened, being a sufficient criterion for the positive evaluation of the Vorlie at least one further secondary significance (1'-n ') occurring within the sub-time interval must be detected, and the opened sub-time interval is preferably maximally 5 minutes, more preferably no more than 1 hour.
[15]
15. The method according to any one of claims 1 to 14, characterized in that as a sufficient condition for the positive evaluation (the presence of a relevant HFl -indicated influencing) the detection of a defined number of secondary significances (1'-n ') is required, preferably at least three, more preferably at least five, per reference time interval, which is preferably between 0.5 and 12 hours, more preferably between 5 and 30 minutes.
[16]
16. Method (1) according to any one of claims 1 and 3 to 15, characterized in that for the synchronization of the HFI and the HRV-FP curve for the considered reference time interval OYn) in each case an equal number of the two courses constituting HFI and HRV -FP Mess¬oder calculation values is generated, preferably by interpolation and / or extrapolationund / or selection of detected HFI and HRV-FP measurement / calculation values, and temporally korrespondonierende HFI and HRV-FP measurement / calculation values each with the same time indices In a preferred embodiment, the reference time interval (T, .n) for both (HFI and HRV-FP) courses is subdivided into a multiplicity of sub-time intervals of a maximum of 30 seconds, particularly preferably of a maximum of 10 seconds , and in each of these sub-time intervals, an HFI measurement / calculation value and an HRV-FP measurement / calculation value are respectively selected, each of which is used as a detection criterion for prim and secondary signatures are most significant, preferably minimum, maximum, median, average or delta values, and the HFI and HRV-FP measurement / calculation values selected in the respective sub-time interval in the sequence checked for correlation.
[17]
17. The method according to any one of claims 1 and 3 to 16, characterized in that for the positive evaluation (the presence of a relevant HFl-indicated influence) at least one, preferably at least three coincident Significanzen (1-n) must be detected per hour or that in the analyzed reference Time interval more than 20%, preferably more than 30%, more preferably more than 50% of the detected in the HFI course primary significances (1-n) with coincident secondary significances (1'-n ') in the HRV-FP curve correlate.
[18]
18. The method according to any one of claims 1 and 3 to 17, characterized in that one or more of the following characteristics are detected in a preferably automated way as the primary significances (1-n) in the HFI curve: - Maximum values or power flow peaks in the case of a graphical representation, as approximately needle-shaped or jagged amplitude deflections, approximately step-shaped course sections, in this case preferably the beginning and / or end sections of an approximately horizontal or a constant power flux density indicative course section of the staircase form and / or the same approximately horizon ¬talen course section itself, - several, preferably more than three, preferably within a time interval of a maximum of 10, 30 or 60 minutes consecutive increases and decreases in the course, - exceeded a power flux density of more than 0.1 mW / m2, preferably more than0, 05 m W / m 2, particularly preferably more than 0.01 mW / m 2, - approximately constant power flux densities or, in the case of a graphical representation, substantially parallel to the time axis extending sections of the HFI curve, - Reductions and / or increases in the HFI curve or the power flux density, which preferably have a difference of more than 10%, particularly preferably more than 30%, compared to immediately preceding and / or following HFI progress sections, where as immediately preceding a period of several, eg 1-30 minutes, preferably from 0-60 seconds, more preferably from 0-10 seconds can be defined.
[19]
Method according to one of Claims 1 and 3 to 18, characterized in that the interval duration, i. the time intervals of successive primary significances (1-n) in the HFI course, and preferably also their intensity, i. subjecting the power density density difference between immediately preceding and / or subsequent HFI history sections to frequency analysis generating statistical characteristics for each of which thresholds or maxima and / or minima are determined, wherein exceeding or undershooting these thresholds is recognized as the primary meta-significance; the occurrence of such a primary meta-significance is checked for a temporal correlation with secondary significances detected in the HRV-FP course, wherein the presence of such a correlation is evaluated as an indication of the presence of a human biologically relevant influencing of the vegetative regulatory mechanisms of the test subject and being used as a statistical characteristic value Preferably, the SDNN (standard deviation of normal-to-normal intervals) related to the interval duration of successive primary significances (1-n) is used.
[20]
A method (1) according to any one of claims 1 and 3 to 19, characterized in that a targeted initiation of a radio frequency immission applied to the test subject is made by one or more dedicated test radio frequency sources and one or more times determined primary significances (1-n) are generated, wherein currently determined HRV-FP curve of the test subject is checked for the primary significance (s) (1-n) coincident secondary significances (1'-n ') and where the initiation of the high-frequency immission preferably in a (on external, ie not caused by the test high-frequency source high frequency immissions related) HFI-course area with approximately constant exposure altitude over a defined period of time or in the presence of external high-frequency Immissions, which are below a defined threshold, vorzugswei¬se below a power flux density of 0.4 mW / m 2, takes place.
[21]
21. Method (1) according to any one of claims 1 and 3 to 20, characterized in that the test subject during one or more time intervals by a preferably surrounding on all sides of the test subject shielding device such. is shielded by electrically conductive metal fabric from external high-frequency immissions of the environment and the Abschirmvorrich¬ device opened either for a defined period of time and thus the test subject again exposed to the external high-frequency immissions or that with the shielding in the time interval of the shield at least one within the shielding is arranged one or more times for a defined period of time, wherein the simultaneously determined HRV-FP profile of the test subject on the / by the test high-frequency source generated primary significance (s) (1 -n) coincident secondary Signifikan¬zen (1 '-n') or on occurring during the period of the opening of the shielding second secondary significances (Tn ') is checked.
[22]
22. Method according to one of claims 1 and 3 to 21, characterized in that from the high-frequency immissions acting on the test subject, a plurality of HFI curves are determined in different frequency band ranges, whereby such frequency-selective HFI curves are in each case detected with at least one HRV detected simultaneously. FP curve and coincident signatures (1-n, Tn ') are detected.
[23]
23. The method according to any one of claims 1 and 3 to 22, characterized in that the HFI course and the HRV-FP course constituent measurement or calculation values in each case to a gradient sequence (GHrv-fp and GHr) are converted, which from a plurality temporally successive Gradients, ie Difference formations between in each case two chronologically successive HFI or HRV-FP measurement or computation values are combined, thus in the case of a graphic representation is shown as a bar graph, the gradient sequences (Ghrv-fp and Ghfi) are analyzed according to the invention in terms of their correlation.
[24]
24. The method according to claim 23, characterized in that the gradient sequences (Ghrv-fp and GHfi) corresponding graphical profiles are created, wherein the HRV-FP Ver¬lauf corresponding gradient sequence (GHrv-Fp) and the HFI gradient corresponding gradient sequence ( Ghfi) are superimposed or juxtaposed, hereby preferably mirrored to a mirror axis parallel or congruent with the time axis, and wherein at least one of the gradient sequences (GHrv-fp and GHfi) can be scaled in such a way that orthogonal to the time axis average elevations of the Gradients are approximately the same size in both graphs.
[25]
25. The method according to claim 24, characterized in that envelopes are applied to the graphic gradients of the gradient sequences (GHrv-fp, Ghfi) to be compared.
[26]
26. The method according to any one of claims 23 to 25, characterized in that between the HRV-FP course corresponding gradient sequence (GHrv-fp) and the HFI gradient corresponding gradient sequence (GHfi) a preferably in percent specified Korrati¬onswert is determined.
[27]
27. The method according to any one of claims 1 to 26, characterized in that as a verification of the analysis, a comparison of the analyzed HRV-FP course with at least one, preferably with at least two of the recited in claims 3 to 11 HRV-FP courses in view of Presence of substantially simultaneous or causally correlated secondary significances (Tn '), or that in at least two, preferably at least three, of the HRV-FP courses or analysis types listed in claims 3 to 11 each have a defined number of secondary significances ( 1'-n ') must be detected.
[28]
Method according to any one of claims 1 to 27, characterized in that, depending on the determined frequency and / or standard deviation (size) of the HFI-induced secondary significances (1'-n '), a gradual, e.g. in percent, the likelihood of the presence of incompatibility of radio frequency emissions is assessed; and / or a degree, e.g. percentile evaluable of the subject's vegetative regulatory ability against HFI-induced secondary significances (1'-n ') is performed on a comparison with HRV-FP comparison data, this comparison data being actual or past measurements on the test subject itself or empirical values of externalHRV-FP test series based and / or a gradual, eg Percentable evaluation based on the speed recorded in the HRV-FP course with which the organism of the test subject after detection of HFI-induced secondary significances (Tn ') in the HRV-FP course returns to a standardized or individually calculated HRV-FP control course, ie to a defined target Values is returned.
[29]
29. Arrangement for the analysis of the interaction of radio frequency electromagnetic emissions (HFI) with vegetative regulatory mechanisms of a test subject according to claim 1, comprising a measuring device for measuring the intensity of high frequency immissions (HFI), preferably determined as power flux density, and an electrocardiogram (EKG) measuring device including measuring electrodes Determining the temporal variation of consecutive heart rate intervals as well as the course of determined, associated with the heart rate variability (HRV) vegetative function parameters (HRV-FP) of the test subject and a processor-controlled evaluation together with this data-connected storage device, on which one or more evaluation algorithms for performing an analysis method according to one of claims 3 to 28, wherein preferably the ECG measuring device and / or the evaluation device and / or the HFI measuring device in one of the test subj ekt applizierbarenGerät are combined or the evaluation is positioned at an external location and by means of the ECG and HFI measuring devices determined records to this external, voresweisweiserverbased evaluation for the implementation of the analysis method according to the invention are überttelbar.
[30]
30. Arrangement for analyzing the interaction of radio frequency electromagnetic emissions (HFI) with vegetative regulation mechanisms of a test subject according to claim 2, comprising an electrocardiogram (EKG) measuring device including measuring electrodes for determining the temporal variation of successive heart rate intervals and the course of said rate variability (HRV) determined therefrom. associated vegetative function parameter (HRV-FP) of the test subject as well as a processor-controlled evaluation device together with this in data-connected storage device on which characteristic, HFl-induced secondary Signifi¬kanzen or history characteristics are stored in the form of digital or analog data and one or more evaluation algorithms for Carrying out an analysis method according to one of claims 3 to 28, wherein preferably the ECG measuring device and / or the evaluation device are combined in a device that can be applied to the test subject or the evaluation device is positioned at an external location and transmits data sets determined by the EKG and HFI measuring devices to this external, preferably server-based evaluation device for carrying out the analysis method according to the invention or for comparison with the characteristic, HFI-induced secondary significances are.

类似技术:

---

公开号 | 公开日 | 专利标题

---

Foster et al.2015|Intrinsic and task-dependent coupling of neuronal population activity in human parietal cortex

---

DE102012218057A1|2014-04-03|DEVICE AND METHOD FOR INVESTIGATING A NARROW INTERACTION BETWEEN DIFFERENT BRAIN SIZES

---

WO2014072461A1|2014-05-15|Method and device for determining vital parameters

---

CN103584872A|2014-02-19|Psychological stress assessment method based on multi-physiological-parameter integration

---

EP1596938A1|2005-11-23|Assembly and method for carrying out magnetotherapy

---

DE102007010353A1|2008-03-13|Diagnosis of sleep apnea

---

DE102016011700A1|2018-03-29|Monitoring of biosignals, in particular electrocardiograms

---

AT516204B1|2016-07-15|Method and arrangement for analyzing the interaction of high frequency electromagnetic emissions with vegetative regulatory mechanisms of a test subject

---

DE112016001171T5|2017-11-30|CVHR-shape measurement apparatus

---

Xu et al.2019|Fatigue EEG feature extraction based on tasks with different physiological states for ubiquitous edge computing

---

EP2799110A1|2014-11-05|Method for operating a bioresonance device

---

EP1106205B1|2008-03-26|Implanted cardiac device with detection of day/night periods

---

DE102015116044A1|2017-03-23|Method and device for quantifying a respiratory sinus arrhythmia and use of such a method or device [

---

DE102009002134B4|2014-10-23|A method for detecting a change in the regulatory system of a human test subject based on its heart rate variability due to changes in the inhomogeneity of a magnetic field applied to the test subject

---

DE102012206293A1|2013-10-17|Method for operating pulse oximeter, involves identifying motion artifacts in measuring signal, where measuring signal is evaluated with respect to arterial oxygen saturation in time intervals of measuring signal

---

DE102014207730A1|2015-10-29|A method, apparatus and computer program product for determining an isoelectric point

---

DE102012112689A1|2013-06-20|Method for determining and promoting health of person by using smartphone, involves performing evaluation by comparing indications with each other, and optically and/or acoustically displaying evaluation result on display unit

---

AT523903A1|2021-12-15|Environmental/health monitoring system for recording electromagnetic fields and influencing vegetative regulation mechanisms of the human organism

---

Kasturiwale et al.2017|Influence of physiological factors on heart rate variability in predicting cardiovascular diseases

---

DE102016015631B4|2019-01-17|Device, system and method for generating biofeedback

---

DE102018107120A1|2018-11-22|Selective activation of parasympathetic nerve in a vagus nerve stimulation device

---

Schumann2011|Fluctuations and synchronization in complex physiological systems

---

CN110192876A|2019-09-03|Based on the lie detecting method for more leading EEG signals kurtosis

---

DE102005029562A1|2006-01-19|Method and device for monitoring cardiac activity

---

DE102011117573A1|2013-05-02|Method and device for the non-invasive monitoring of pain intensities and the determination of depths of anesthesia during medical interventions

同族专利:

---

公开号 | 公开日

---

DE102015011213A1|2016-03-03|

---

AT516204B1|2016-07-15|

引用文献:

---

公开号 | 申请日 | 公开日 | 申请人 | 专利标题

---

WO2003096893A1|2002-05-17|2003-11-27|Iembio Co., Ltd.|Portable heart rate variability based health monitoring system having electromagnetic field sensor built in|

---

RU2303392C1|2005-12-26|2007-07-27|Ярославский государственный технический университет|Method of estimation of degree of influence of electromagnet fields onto human organism|

---

WO2010112503A1|2009-04-02|2010-10-07|Elisabeth Plank|Use of the heart rate variability change to correlate magnetic field changes with physiological sensitivity and method therefor|

---

JP3946108B2|2002-08-27|2007-07-18|パイオニア株式会社|Heart rate variability analyzer, heart rate variability analysis method, and heart rate variability analysis program|

---

DE102006039957B4|2006-08-25|2012-08-16|Biosign Gmbh|Method for evaluating heart rate variability|CN106326644B|2016-08-16|2019-05-17|沈阳东软熙康医疗系统有限公司|A kind of computing device of heart rate variability parameter and fatigue strength index|

---

AT523903A1|2020-06-12|2021-12-15|Matissek Msc Michael|Environmental/health monitoring system for recording electromagnetic fields and influencing vegetative regulation mechanisms of the human organism|

法律状态:

优先权:

---

申请号 | 申请日 | 专利标题

---

ATA666/2014A|AT516204B1|2014-08-29|2014-08-29|Method and arrangement for analyzing the interaction of high frequency electromagnetic emissions with vegetative regulatory mechanisms of a test subject|ATA666/2014A| AT516204B1|2014-08-29|2014-08-29|Method and arrangement for analyzing the interaction of high frequency electromagnetic emissions with vegetative regulatory mechanisms of a test subject|

---

DE102015011213.9A| DE102015011213A1|2014-08-29|2015-08-28|Method for analyzing the interaction of high frequency electromagnetic emissions with vegetative regulatory mechanisms of a test subject|